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1. SUMMARY 


The cluster approach has been applied to the study of the 
electronic structure of cathode surfaces and high-temperature 
superconductors. The reported calculations are based on the 
relativistic Xa scattered-wave method of computational quantum 
chemistry. The study of surface-dipole properties and dipole- 
dipole interactions for low- work -function surfaces has enabled, for 
the first time, a systematic theoretical investigation of work- 
function/coverage curves and a determination of minimum work 
functions for different model cathode surfaces. 

Surface-dipole properties have been obtained for two series of 
calculations. In the first series, BaO adsorbed on monolayer films 
of Re, Os, Ir, and Pt on otherwise unchanged bcc W(100) has been 
investigated to study electronegativity and charge -transfer effects 
for different coating metals. In the second series, BaO on bcc 
W(110) , hep Os (1010), and fee Pt(lll) has been investigated to 
study the effects of different substrate crystal structures. It 
was found that the surface-dipole properties behave monotonically 
for the monolayer films, while a minimum in the surface dipole_and 
a maximum in the dipole energy is obtained for BaO/Os (1010) . 
Alloying with tungsten increases_both the surface dipole and barium 
binding energy for the hep Os (1010) substrate, while it has a small 
effect in the case of the cubic substrates. 

Extended cluster calculations have been carried out for the 
study of adsorbate interactions and the dependence of the dipole 
properties on the BaO surface coverage. Larger clusters are found 
to provide systematically improved results for the surface 
electronic structure and the adsorbate binding energy. At higher 
coverages, the calculated binding energies for Ba on B-type cathode 
surfaces of 4. 6-5. 4 eV compare very well with the observed energy 
of desorption of 4.8 eV. The calculated binding energies for Ba 
and Cs on W(100) of 3.6 and 3.0 eV, respectively, are in excellent 
agreement with the experimental energies of desorption of 3.8 and 
2.8 eV. 


:K. : The depolarization model of interacting surface dipoles has 
been employed to determine, from first principles, minimum work 
functions for a variety of low-work- function surfaces relative to 
BaO/W. The two important microscopic properties, the initial 
surface dipoles at low coverages and their polarizabilities, or 
depolarization, have been determined from calculations for clusters 
with individual and interacting surface dipoles. Consistent and 
qualitatively accurate results have been obtained for minimum work 
functions_ relative to BaO/W (2.0 eV) , namely 1.8 eV for 
BaO/Os (1010) , 2 .4 eV for BaO/Pt (111) , 2 . 6 eV for Ba/W (100) , and 1 . 4 
eV for Cs/W (100) . 
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Comparisons with theoretical results have been made for the 
Ba, 0, and W electronic energies observed in photoelectron spectra 
of Ba0/W(100), the oxygen- induced features of 0/W(100) in inverse 
photoemission spectra, and spectroscopic constants of BaO. 
Different orientations of BaO on the surface have been considered. 
For the Ba-0 interatomic distance on W, single-adsorbate cluster 
calculations favor the value of 2.3 A suggested by Shih et al . 
(1988) , but the question of the orientation of BaO on operating 
cathodes surfaces remains unresolved . 

The electronic structure of the conduction and valence bands 
and the electronic core levels of the 90K superconductor YBa 2 Cu 3 0 7 
have been studied with fully relativistic calculations. Different 
clusters including between 1 and 4 Cu-0 layers and up to 39 atoms 
have been considered. The generated densities of states for 
clusters of the size of a unit cell are found to be in good 
agreement with solid-state calculations and experimental results. 
Projected densities of states show that the characteristic double 
maximum observed in photoemission spectra just below the Fermi 
level is composed of Cu states at its center and 0 states more to 
its sides. The high oxygen density of states near the Fermi level 
and the calculated charge distribution point to substantial 
vacancies in the 0 2p shell. Large covalent contributions in the 
chemical bonds are found which are expected to be characteristic 
for the surface of the material. 

Core -electron binding energies have been calculated for all 
core states of the two inequivalent Cu atoms and the four 
inequivalent 0 atoms, as well as the 3d states of Y and Ba. Final- 
state effects in the photoemission process and spin-orbit 
interactions in the valence band have been investigated for the 
first time. While final-state effects are essential for the 
calculation of accurate core-level energies, they are found to be 
small (0.2 eV) for the localized Cu 3d states in the valence band 
and do not significantly affect its location relative to the Fermi 
level. The Cu 3d spin- orbit splittings are 0.3 eV. 

The results for BaO/W cathode surfaces and the high- 
temperature superconductor YBa 2 Cu 3 0 7 have been presented at several 
scientific conferences and published in the IEEE Transactions on 
Electron Devices and in Physica C (Muller, 1989) . 
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2 . INTRODUCTION 


It is well known that the electron emission from impregnated 
barium dispenser cathodes can significantly be improved by coating 
the tungsten surface with a platinum-group transition metal, like 
Re, Os, Os-Ru, or Ir. Similar effects can be achieved by using 
mixed-metal matrices instead of a pure tungsten matrix. During the 
activation and early life of the cathode, interdiffusion of W and 
the other metal (s) occurs and an alloy surface is formed, with a 
monolayer of adsorbed Ba and 0 from the impregnant . The best 
coating metal is Os, followed by Ru, Ir, and Re, while, e.g., Pt 
gives rise to a degradation of the emission properties. 

The improved emission is caused by a lowering of the work 
function from about 2 . 0 eV for tungsten (B-type) cathodes to 1.8 eV 
for alloy (M-type or mixed-metal) cathodes. Pt coating, on the 
other hand, increases the work function to 2.2 eV (Norman et al . , 
1987) . A clean W surface, for comparison, has a work function of 
4.5 eV. The decrease from 4.5 eV for the clean substrate to around 
2 eV for cathode surfaces is caused by the adsorption of Ba and 0 
and attributed to the formation of a surface dipole that lowers the 
electrostatic potential for electron emission. The effect due to 
the alloy formation is comparatively small, less than 10 percent of 
the total effect, but nevertheless technologically very important 
because a work function variation of 0.2 eV changes the emission- 
current density (at the same temperature) by one order of 
magnitude . 

Several attempts have been made to explain the variation of 
the work function and emission that is observed for cathodes with 
different coating materials. From alkali-metal adsorption studies 
it is known that usually the substrates with the highest work 
functions give rise to the lowest work functions after adsorption. 
By analogy, Zalm and van Stratum (1966) explained the emission 
enhancement for Re, Os, Ru, and Ir by the higher work functions of 
these metals as compared to W. One would therefore expect that a 
coating with platinum, which has an even higher work function, 
should lead to a further improvement. This is, however, not the 
case . 


Norman et al . (1987) have investigated the bonding of Ba and 
O on W, W-Os, and W-Pt cathodes, using surface extended x-ray- 
absorption fine-structure (SEXAFS) experiments. They find that Ba 
is coordinated with one O atom for BaO/W and BaO/Pt-W, but coordi- 
nated with two 0 atoms for BaO/W-Os, and tentatively explain the 
lowering of the work function by the larger charge transfer from Ba 
being bridge-bonded to two O atoms. SEXAFS experiments performed 
by Shih et al . (1988) for the B-type surface, however, contradict 
the assignments of the earlier experiment, and they propose a 
bridge-bonded structure for B-type surfaces. Calculations by 
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Hemstreet et al . (1989) for Ba0/W(100) appear to confirm Shih' s 
results, but the calculated work functions are significantly higher 
than expected. 

The purpose of the present theoretical study is to 
investigate, with first-principles cluster calculations, possible 
reasons for the observed emission characteristics of different low- 
work- function surfaces, and to calculate the work functions of 
different materials. 

This report is organized as follows: The results for models 
of B- and M-type cathode surfaces are presented in Section 3. The 
electronic structure of 0/W and ab initio results for BaO are 
described in Sections 4 and 5, respectively. An application of the 
relativistic cluster approach to the recently discovered high- 
temperature superconductor YBa 2 Cu 3 0 7 is presented in Section 6 . The 
computational methods are briefly discussed in Section 7, and the 
conclusions are summarized in Section 8. 


3. CATHODE SURFACES 


In this section, the surface-dipole properties, work 
functions, and electronic structure of low-work- function surfaces 
are described. The purpose of these investigations is to identify 
the microscopic materials parameters that determine the different 
work functions and significant variations in emission enhancement 
for B- and M-type cathode surfaces. 


3.1 Surface-Dipole Properties 


The surface -dipole properties are compared for two series of 
calculations of model cathode surfaces. The first series involves 
BaO adsorbed on monolayer films of Re, Os, Ir, and Pt on an 
otherwise unchanged bcc W(100) substrate. In the second series, 
effects of the different crystal structures of the coating 
materials and alloys are investigated. In this series, the 
dominant crystal faces as present on actual cathode surfaces are 
employed . 


3.1.1 Monolayer Films 

In order to investigate the differences in the electronic 
effects of W, Re, Os, Ir, and Pt without the complications arising 
from varying crystal structures, different crystal planes, and 
uncertain adsorption geometries, monolayer-film properties are 
discussed first. These results reflect the changes in the 
electronic structure and surface-dipole properties induced by the 
presence of the coating metal in otherwise identical geometrical 
environments . 

The substrate cluster used is the unit cell (a W 9 cube) of 
body-centered-cubic tungsten with Ba adsorbed on top of 0 at the 
fourfold-hollow site of the (100) face. For the monolayer films, 
the W surface layer is replaced with Re, Os, Ir, or Pt . Because 
there are no adsorbate interactions present in these clusters, the 
results are to be considered as representative of low-coverage 
conditions . 

The calculated properties include the charge distribution with 
net atomic charges q-t at locations r^ the resulting 


surface dipole 

M = £ i , 

(1) 

dipole length 

R di P = p/q(A) , 

(2) 

and dipole energy 

E dip = q 2 (A)/R dip = q 3 (A) / [i . 
5 

(3) 



q (A) is the net charge on the adsorbate (typically Ba) . E dip 
is a measure of the attractive contribution to the adsorbate 
binding energy and may be compared with experimental energies of 
desorption. Atomic units are used, i.e., m = e = h/27r =1. In 
these units q is expressed in electrons (e) , r and R in Bohrs (la 0 
= 0 . 5292 A) , n in a.u. (= 2 . 542 Debyes) , and E in Hartrees (= 27 .21 
eV) . 


The results for the clusters BaO/M^W^ with a surface layer of 
M = W, Re, Os, Ir, and Pt are given in Table I. The surface -dipole 
properties were calculated using the quasirelativistic scattered- 
wave Xa approach. They are in good agreement with fully 
relativistic results (see Section 7) . 


TABLE I. 

Dipole properties for BaO/M 4 W 1 W 4 (bcc,100) clusters with M = W, Re, 
Os, Ir, and Pt from quasirelativistic Xa calculations. q (in 
electrons) denotes net atomic charges, n dipole moment, R dip dipole 
length, and E Ba dipole energy. R(Ba-O) = 4.9 a 0 «* 2.6 A and 
R (Ba-surface) = 6.0 a 0 « 3.2 A (Tuck and Norman et al . , 1986 and 
1987) . 


Property M 

W 

Re 

Os 

Ir 

Pt 

q (Ba) 

+ 1.50 

+ 1.53 

+1.57 

+ 1.58 

+ 1.52 

q (0) 

-0.15 

-0.17 

-0.15 

-0.10 

-0.25 

q (Substrate) 

-1.35 

-1.36 

-1,42 

-1.48 

-1.27 

Electronegativity 

1.40 

1.46 

1.52 

1.55 

1.44 

q(M x ) 

-0.08 

-0.17 

-0.25 

-0.34 

-0.27 

q(W n ) 

-0.00 

+ 0.06 

+ 0.12 

+ 0.24 

+ 0.18 

q(w IIX ) 

-0.25 

-0.19 

-0.13 

-0 . 09 

-0.09 

H (Debyes) 

38.4 

35.4 

31.1 

27.3 

25 . 8 

Rdip (A) 

5.3 

4.8 

4.1 

3.6 

3.5 

E Ba (eV) 

6.1 

7.0 

8.6 

10.0 

9.4 


Barium is found to be heavily oxidized with a net charge of 
1.5-1. 6 electrons. Only about 10-20% of this charge has been 
transferred to oxygen, while the major fraction has been 
transferred to the substrate. A comparison between the different 
clusters shows that Ba has the highest oxidation level for the Ir 
monolayer film. Experimentally, Haas et al . (1983) have found that 
the oxidation level increases for BaO adsorbed on W to W-Os and Ir, 
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whic:h is in excellent agreement with the present results. This 
observed behavior closely parallels and has its origin in the 

metals 1 (cf Tabie the different coating 

The largest charge redistributions are found among the 
substrate atoms. With increasing electronegativity of the atoms in 
the first substrate layer, this layer acquires more and more 
electrons from the underlying W atoms. Because of the strong 
interaction of 0 with the second-layer W atom in this geometry, the 
second layer becomes slightly oxidized. y 

The large charge transfer from Ba to the substrate gives rise 
to very large surface dipole moments. The calculated dipole 
moments decrease monotonically in the series from BaO on W to Pt-W 
in spite of the maximum in the net charge on Ba for the Ir-W 
system, because the charge from Ba is increasingly transferred to 
f lrst su ^ strate layer at the expense of mainly the third layer 
T ?^ S v. le ^ ds to a si 9 nifi cant reduction in the dipole length R ’ 
which, . because of a comparatively stable net charge on Ba 
determines the observed trend in the surface dipoles. ' 

^ dipole energy increases for the monolayer clusters from 

^ • 1 eV f° r W to 10.0 eV for Ir and then falls back to 9.4 eV for 
Pt This trend is expected to reflect the variation of the Ba 
binding (or desorption) energies at very low coverages, and is 
qualitatively in agreement with the observed increase of the Ba 
desorption energy from 4.8 to 5.4 eV for B- to M-type cathodes 
surfaces (Forman, 1984 and 1987) . 

The results, for the monolayer clusters are summarized in Fig 
1, which shows the monotonic decrease of the surface dipoles and 
the increase m the dipole energies. Both trends are governed by 
the increasing electronegativity of the refractory metal (up to Ir) 
in the first substrate layer. The increasing transfer of charqe 
from the underlying W atoms to the first layer reduces the Ba- 
substrate dipole and explains the decrease in the surface dipoles 
(AU • Because the Ba binding energy is proportional to 1 /u, the Ba 
binding energies increase accordingly. The maximum for Ir-W 
results from the fact that Ir has the highest electronegativity of 
the. refractory metals, and as a consequence the net charge (q) on 
Ba is a maximum. Because the binding energy is proportional to the 
second power of q, the effect of the further decrease in R,. (and 

ty. ^? r Pt-W is more than compensated for and a maximum in the 
binding energy is obtained. 



FIGURE 1. 



BaOonW Re/W Os/W Ir/W Pt/W 

(bcc) (bcc) (bcc) (bcc) (bcc) 
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3.1.2 Effects of Crystal Structure 

For the study of the effects of different crystal structures 
the dominant faces as present on actual cathode surfaces are 
employed. These are the (110) face of bcc W, the (1010) face of 
hep Os-W, and the (111) face of fee Pt-W. Lattice constants of 
3.16 A for W, 2.74 and 4.32 A for Os, and 3.92 A for Pt are used. 

Ba is considered to be adsorbed on top of O at a hollow site 
in most cases. Despite the recent verification Tshih et al . 1991) 
that Ba and 0 are adsorbed at room temperature at different sites 
(and not on top) for BaO/W(100) with a Ba-0 distance of 3.2 A, the 
distance on actual cathode surfaces is so much shorter ( 2 . 3 - 2. 6 A) 
that on-top adsorption is still the most reasonable assumption. 
For the Ba-W(Os) distances the values suggested by Tuck and Norman 
et al . (1986 and 1987) were used. The different clusters and Ba-O 
geometries that were employed for W, Os, Pt, and alloyed substrates 
are described in Table II. The clusters are shown in Fig. 2. 


TABLE II. 


Substrates, crystal faces, clusters, Ba-0 orientation, and 
interatomic distances used for the study of cathode surfaces (1 a„ 
= 0.5292 A) . 


Substrate 

Face 

Cluster 

Ba-0 

R (Ba-O) 

R (Ba-surface) 

bcc W(-Os) 

( 110 ) 

BaO/ (W-Os ) 9 

on- top 

o 

m 

Lfi 

Ch 

7.39 

a o 

hep Os ( -W) 

( 1010 ) 

BaO/ (Os-W ) 9 

on- top 

4.90 a 0 

6 . 90 

a o 

hep Os ( -W) 

( 1010 ) 

Ba0 2 / (Os-W) 9 

bridge 

4.90 a 0 

5 .20 

a o 

fee Pt(-W) 

( 111 ) 

BaO/ (Pt-W ) 7 

on- top 

5.10 a 0 

6 . 90 

a o 


In comparison to the monolayer- film results, the trends 
obtained for the different substrate crystal structures are 
completely different. As seen in Table III and shown in Fig. 3, 
the dipole moment exhibits a minimum and the dipole energy a 
maximum for the Os hep substrate. The maximum in the dipole energy 
is consistent with the observed increase in the Ba binding energy 
from B- to M-type cathode surfaces, but the variation in the dipole 
moment is completely opposite to expectation. A large (initial) 
surface dipole is usually indicative of a lower effective work 
function. The calculated dipole for Os (hep) , however, is smaller 
and the dipole for Pt (fee) is larger than that for W (bcc) , while 
an Os surface is a good emitter and a Pt surface a poor emitter. 
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FIGURE 2. 


Clusters for the study of BaO on W(110), Os (1010), and Pt(lll) . 



BOQ/WgCLLO) 

bcc 


BaO/OSgdDlO) 

hep 


BaO/Pt 7 (lH) 

fee 


TABLE III. 

Dipole properties for BaO/W 9 (bcc, 110) , BaO/Os 9 (hep, 1010) , and 
BaO/Pt 7 (fee, 111) clusters from quasirelativistic Xa calculations, 
q (in electrons) denotes net atomic charges, n dipole moment, R dip 
dipole length, and E Ba dipole energy. 




Property 

BaO/W 9 

BaO/Os 9 

BaO/Pt 7 

q (Ba) 

+1.48 

+1.57 

+ 1.48 

q(O) 

-0.49 

-0.42 

-0 . 12 

q (Substrate) 

-0.99 

-1.15 

-1.36 

q(M x ) 

+ 0.00 

-0.22 

-0 . 04 

q(M„) 

-0.66 

+0.17 

-0.67 

q(Mm) 

-0.08 

-0.01 

-0.23 

q (Mj V ) 


-0.21 


/x (Debyes) 

39.6 

33.3 

45.4 

Rdip (A) 

5.6 

4 . 4 

6.4 

Esa (eV) 

5.7 

8.0 

4.9 
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Dipole Energy (eV) Dipole Moment (Debyes) 


FIGURE 3. 


Surface dipoles and Ba binding energies for BaO on W(bcc 110) 
Os (hep, 1010) , and Pt (fee, 111) [full symbols] 
and their alloys [open symbols] . 


SURFACE DIPOLES 



Ba BINDING ENERGIES 



ii 






In order to verify these surprising results for more realistic 
substrates, quasirelativist ic Xa calculations have been carried out 
for randomly alloyed substrate clusters for each crystal structure. 
Tables IV- VI I show that the charge distributions and dipole 
properties for the alloys with the cubic crystal structure are 
similar to the results for the pure substrates. The largest 
effects due to alloying are found for the Os hep crystal structure, 
where the net charge on Ba increases by about 0.2 electrons and the 
net charge on 0 doubles in most cases. The surface dipole for the 
Os-W hep alloy is almost as large as for the W bcc substrate, due 
to the larger charge transfer but smaller dipole length. 

The net charge on Ba and the surface dipole decreases signifi- 
cantly if Ba is bridge-bonded to two O atoms and closer to the 
surface (cf . Tables V and VII) . Alloying again increases the 
dipole moment and, to a larger extent, the dipole energy. The 
bridge -bonded results are shown in Fig. 3 by squares. 

The effects arising from the different substrate crystal 
structures of W (bcc), Os (hep), and Pt (fee) and their alloys, as 
compared to the monolayer films, can be explained as follows: 

In Table III it is seen that the net charge on Ba runs through 
a maximum for Os, which is very similar to the monolayer- film 
results. The dipole moment, however, shows a minimum for the pure 
Os substrate, while the dipole moment behaved monotonically for the 
monolayer films. The minimum in the dipole for Os (in spite of the 
maximum in the charge transfer from Ba) has its origin in the 
different crystal structures for W, Os, and Pt. For the present 
faces, the interlayer spacing is smallest for Os and largest for 
Pt. This obviously has the effect of a better screening of the 
positive charge outside the surface, with the consequence that the 
negative charge balancing the positive charge from Ba is located 
closer to the surface in the case of Os, while more image-like 
charges are present in the case of W and particularly Pt . An 
increase in the negative charge on 0 is also observed for the W and 
Os substrates relative to the monolayer-film results. The dipole 
energies behave opposite to the dipole moments as before, and both 
trends are again governed by the length of the surface dipole, 
which is affected by the spacing and resultant screening in the 
substrate layers . 

In conclusion, while it is found that the dipole moments and 
energies vary monotonically for the monolayer-film clusters, the 
dipole moment exhibits a minimum and the dipole energy a maximum 
for the hep crystal structure of the Os substrate, relative to bcc 
W and fee Pt . Alloying with W increases both the dipole moment and 
energy for the hep substrate . 
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TABLE IV. 


Dipole properties for BaO/(W-Os) 9 (bcc,110) clusters. q (in 
electrons) denotes net atomic charges, /x dipole moment, R dip dipole 
length, and E Ba dipole energy. 


Property 

BaO/W 9 

BaO/W 5 Os 4 

BaO/W 4 Os 5 

q (Ba) 

+ 1.48 

+ 1.47 

+ 1.48 

q (0) 

-0.49 

-0.47 

-0.46 

q (Substrate) 

-0 . 99 

-1.00 

-1.02 

q(W I ,0s I ) 

+ 0.00 

-0.05,-0.04 

-0 . 01, +0 . 02 

q(W IIf 0s„) 

-0.66 

-0.54, 

, -0.78 

<3 ^Hi/ OSm) 

-0 . 08 

+0.20,-0.34 

+0.23, -0.36 

/x (Debyes) 

39.6 

38.0 

39.8 

R dip (A) 

5.6 

5.4 

5.6 

E Ba (eV) 

5.7 

5.8 

5.6 


TABLE V. 

Dipole properties for BaO/(Os-W) 9 (hep, 1010) 
electrons) denotes net atomic charges, /x dipole 
length, and E Ba dipole energy. 

clusters. q (in 

moment, R dip dipole 

Property 

BaO/Os 9 

BaO/ Os 5 W 4 

BaO/Os 4 W 5 

q (Ba) 

+ 1 . 57 

+1.77 

+ 1.75 

q (0) 

-0.42 

-0.98 

-0 . 97 

q (Substrate) 

-1.15 

- 0.79 

-0.78 

q (Os I7 Wj) 

-0.22 

-0.19, +0.14 

-0.15, +0.18 

q (OSjj,Wjj) 

+ 0 . 17 

, +0.20 

-0.19, 

q (OSjjj, W riI ) 

-0 . 01 

-0.12, +0.23 

-0.22, +0.12 

q (OSjy , W IV ) 

-0.21 

-0.49, 

, -0.27 

/x (Debyes) 

33.3 

38.9 

35 . 9 

Rdip (A) 

4.4 

4.6 

4.3 

E Ba (eV) 

8 . 0 

9.9 

10.3 






TABLE VI . 


Dipole properties for BaO/(Pt-W) 7 (fee, 111) clusters. q (in 
electrons) denotes net atomic charges, fi dipole moment, R dip dipole 
length, and E Ba dipole energy. 


Property 

BaO/Pt 7 

BaO/Pt 4 W 3 

BaO/Pt 3 W 4 

q (Ba) 

+ 1.48 

+ 1.56 

+ 1.48 

q (0) 

-0 . 12 

-0.22 

-0 . 06 

q (Substrate) 

-1.36 

-1.34 

-1.42 

<3 ( Pt i / Wj ) 

-0 . 04 

1 

o 

o 

* 

+ 

o 

o 

VO 

-0.20, +0.00 

q ( Pt j X , Wjj) 

-0.67 

-0.76 

-0.63, 

<3 (Ptiu , w m ) 

-0.23 

-0.22,-0.05 

-0.43, -0.09 

fi (Debyes) 

45.4 

45.1 

45 . 6 

R dip (A) 

6.4 

6.0 

6.4 

E Ba (eV) 

4 . 9 

5.8 

4 . 9 


TABLE VII. 

Dipole properties for bridge-bonded Ba0 2 /(0s-W) 9 (hcp,1010) 
clusters. q (in electrons) denotes net atomic charges, fi i dipole 
moment, R dip dipole length, and E Ba dipole energy. 

Property 

Ba0 2 /0s 9 


Ba0 2 /0s 5 W 4 

Ba0 2 /0s 4 W 5 

q (Ba) 

+ 1.22 


+ 1.66 

+ 1.67 

q (O) 

-0.41 


-0 . 90 

-0.44 

q (Substrate) 

-0 . 81 


-0.76 

-1.23 

qfOs^W!) 

+ 0.07 


+0.11, +0.42 

-0.05, +0.31 

q(0s II ,W II ) 

-0.11 


,-0.16 

-0.61, ---- 

q(os m ,w in ) 

-0.08 


-0.16, +0.22 

-0.26, +0.08 

q (Os IV/ W IV ) 

-0.20 


-0.41,- 

, -0.26 

H (Debyes) 

22.6 


27.3 

30.7 

Rdip (A) 

3 . 9 


3.4 

3 . 8 

E Ba (eV) 

5.6 


11.6 

10.5 
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3.1.3 Dipole Interactions 


Two different types of dipole interactions are investigated: 
(a) adsorbate interactions on a surface and (b) dipole interactions 
on opposite sites of a symmetrical slab as often used in solid- 
state calculations. 

(a) Adsorbate Interactions 

First, a comparison is made for fully relativistic_ results for 
vertically and bridge-bonded BaO on W(100) and Os (1010) (Muller, 
IEEE, 1989) . Subsequently, the dipole properties are investigated 
as a function of surface coverage with more extended clusters by 
using the quasirelativistic Xa scattered-wave approach. 

The clusters modeling vertical (BaO/W 9 , Ba 2 0 2 /W 14 ) and 
bridge-bonded (Ba0 2 /W 14 and Ba0 2 /0s 9 ) arrangements are shown in Fig. 
4 . Interatomic distances as proposed by Tuck and Norman et al . 
(1986 and 1987) are used, i.e., R(Ba-O) = 2.6 A, R(Ba-W) = 3.9 A, 
and R(Ba-Os) = 3.75 A. For Ba0 2 /W 14 , the positions of the oxygen 
atoms are kept the same as in Ba 2 0 2 /W 14 and Ba is moved to the bridge 
position such that R(Ba-O) = 2.6 A is maintained. 


FIGURE 4. 


Cluster models for vertical and bridge-bonded BaO on W and Os. 



hep 





Ba02/0s 9 
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In Table VIII and Fig. 5 net charges, surface dipoles, and Ba 
binding energies are compared for the different clusters at the 
fully rel a tivi stic level. The inte racti on of the adsorbate -induced 
dipoles in Ba 2 0 2 /W 14 leads relative to BaO/W 9 to a substantial 
reduction in the charge transfer, dipole moment, and Ba binding 
energy per adsorbate. The dipole moment decreases from 37 to 25 
Debyes and the binding energy from 6.4 to 5.1 eV. The latter value 
compares very favorably with the experimental energy of 4.8±0.2 eV 
for Ba desorption from B-type cathode surfaces (Forman, 1984) . 

For bridge -bonded barium in Ba0 2 /W 14 or Ba0 2 /0s 9 the dipole 
moments also decrease strongly, but the Ba binding energies 
increase relative to BaO/W 9 . The Ba overlayer in the bridge-bonded 
geometry could therefore be more stable and would lead to a minimum 
in the work function at a higher coverage than for the vertical 
orientation because of smaller dipole-dipole interactions. The 
vertical orientation, on the other hand, would lead to an initially 
steeper slope and possibly lower work function at a smaller 
coverage. These results are in qualitative agreement with the 
solid-state calculations of Hemstreet et al. (1989) who find that, 
for low temperatures, bridge-bonded BaO is more stable while 
vertically-bonded BaO produces the lower work function. The 
calculated work functions of 2.7 and 2.6 eV are however 
significantly higher than observed experimentally. 

TABLE VIII. 

Ba 6sp and 5d orbital populations p (in electrons) , net atomic 
charges q (in electrons) , dipole moments n, and Ba dipole energies 
E Ba for BaO /W (Os) clusters from fully relativistic Xa calculations. 


Property 

BaO/W 9 

B3>2®2/^14 

Ba0 2 /W 14 

Ba0 2 /0s 9 

p (Ba, 6sp) 

0.15 

0.16*2 

0.13 

0.23 

p (Ba, 5d) 

0.34 

0.62*2 

0.50 

0.46 

q (Ba) 

+ 1.51 

+ 1.22*2 

+1.37 

+ 1.31 

q (0) 

-0.24 

-0.35*2 

-0.33*2 

-0.36*2 

q (M x ) 

+0.06 

+ 0.01 a 

+ 0 . 86 b 

+ 0 . 05 

q(M„) 

-0 . 33 

-0.89 

-0.85 

-0 . 15 

q(M III<IV ) 

-1.00 

-0 . 86 a 

-0 . 72 b 

-0.49 

/x (Debyes) 

37.5 

24 . 6*2 

26.7 

19 . 7 

E Ba (eV) 

6.4 

5.1*2 

6.7 

7 . 9 

Ba Coverage 

1/4 

1/3 

" i/6 

1/4 

a Wjj-0.18 center, +0. 

19 periph. 

, W m :-0.38 

center, -0.48 

periph. 

b Wji-0.08 center, +0. 

94 periph. 

, W m :-0.32 

center, -0.40 

periph . 
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More extensive calculations for the study of adsorbate 
interactions have been carried out for BaO, Ba, and Cs on W(100), 
and BaO on other substrates by using the quasirelat ivistic 
approach. The surface -dipole results are given in Table IX. 

Fig. 6 shows the clusters used for the Ba0/W(10 0) study, which 
contain 2 or 3 adsorbed BaO molecules and up to 24 W substrate 
atoms. The formal cluster surface coverages (BaO : W surface) 
range from 0.20 to 0.38. The calculated dipole moments and dipole 
energies per adsorbed BaO for BaO/W(100) are plotted in Fig. 7 as 
a function of formal coverage. At coverages of 0.25 and less, the 
dipole moment and energy reach a constant value which represents 
initial adsorption conditions. At higher coverages, repulsive 
dipole-dipole interactions give rise to a large depolarization. 
For a formal coverage of 0.38, the calculated dipole energy for 
BaO/W(100) is close to the experimental value for Ba desorption 
from B-type cathode surfaces. 

Similar calculations have been carried out for Ba and Cs on 
W (100) with Ba- and Cs-surface distances of 4.5 and 5.2 a 0 , 
respectively, and for BaO on W(110), Os (1010), and Pt(lll) . As for 
W (100) , the larger clusters for the other substrates are obtained 
by combining two or more single-adsorbate clusters (cf. Figs. 2 & 
4). The _substrate clusters used for W(110) are W 30 , W 23 , and W 16 ; 
for Os (1010) : Os 21 and Os 15 ; and for Pt(lll) : Pt 16 and Pt 12 . 

For W(110) , the W 9 substrate clusters are combined such that 
R(Ba-Ba) = 5.98 a 0 for the smallest combined unit, which is the 


TABLE IX. 

Dipole moments n (Debyes) as a function of surface coverage for 
BaO, Ba, and Cs on W (100) , and BaO on W( 110) , Os (1010), and Pt (111) 
from quasirelativistic Xa calculations. q 0 are the net charges on 
Ba (or Cs) for the clusters with the smallest surface coverage. 


System 

q 0 

Cov = 0.15 

0.20 

0.25 

0.29 

0.33 

0.38 

0.40 

BaO/W (100) 

1.62 

_ 

30 . 0 

31.2 



26.4 

18 . 0 



Ba/W (100) 

1.30 

— 

20.5 

21.8 


15 . 9 

12.8 

— 

Cs/W (100) 

0 . 92 

— 

— 

18.4 

— 

19 . 1 

16.0 

“ 

BaO/W (110) 

1.60 

42.2 

43 .7 

— 

33.7 

25 . 7 

— 

— 

BaO/Os (1010) 

1.71 

— 

— 

40.5 

— 

32.4 

25.7 

— 

Ba0 2 /0s (1010) 

1.26 

— 

— 

16.6 

— 

17.7 

14 . 1 

— 

BaO/Pt (111) 

1.82 

— 

— 

— 

— 

51 . 9 

— 

34.5 
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FIGURE 7. 


Dipole moments and dipole energies for Ba0/W(100) 
as a function of formal surface coverage. 
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FIGURE 8. 


Dipole moments as a 'function of surface coverage 
for BaO on W(bcc), Os (hep) , and Pt(fcc). 



Surface Coverage 
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same as for the W(100) systems. For the vertical orientation of 
BaO on Os, the 0s 9 substrate clusters are combined such that R(Ba- 
Ba) = 5.18 a 0 , and for the bridge-bonded orientation such that 

R(Ba-Ba) = 8.16 a q . For Pt(lll), portions of the Pt 7 cluster are 
added such that R(Ba-Ba) = 6.04 a 0 . 

The results for BaO on W(100), Os(lOlO), and Pt(lll) are 
compared in Fig. 8. The surface dipoles on Pt exhibit the largest 
interaction (or depolarization) , as is reflected in the steeper 
slope of the dipole/coverage curve at .higher, coverages. The 
depolarization of vertically-bonded BaO is similar on Os and W, 
while bridge-bonded Ba0 2 on Os shows the smallest depolarization. 
These data will form the basis for the determination of work 
function/coverage curves in Sec. 3.2 below. 


(b) Symmetrical Slab 

In solid-state band- 
structure calculations of 
surface properties, the missing 
periodicity in the "z" 
direction perpendicular to the 
surface is usually enforced by 
representing the substrate by a 
slab with a thickness of a few 
atomic layers, and adsorbates 
positioned on both sides of the 
thin substrate film. The 
electronic structure of the 
combined overlayer- substrate 
system is then calculated for a 
periodical arrangement in which 
the slab is repeated in the z 
direction at sufficiently large 
separations to prevent 
interactions between 
neighboring slabs . 

The magnitude of the 
artificial interaction of the 
surface dipoles on both sides 
of an individual slab can be 
evaluated with cluster 
calculations by comparing the 
results for, e.g., a BaO/W 9 
surface cluster with three 
substrate layers and the 
symmetrical Ba 2 0 2 /W 14 slab- 
geometry cluster shown in Fig. 
9. 


FIGURE 9. 

Symmetrical slab -geometry 
cluster Ba 2 0 2 /W 14 . 
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Symmetrical slab-geometry clusters have been investigated for 
Cs, Ba, and BaO on W using the fully relativistic Xa approach. The 
data are compared with earlier surface cluster results (Muller, 
1988) in Table X. As a general result, it is found that the net 
charges q on the adsorbate atoms decrease in all cases. The charge 
on the third, central W layer also decreases significantly, while 
the charge on the second W layer stays almost constant and the 
charge on the first W layer increases. The net charge on the 
adsorbate atoms decreases by about 20% for Cs and 10% for Ba . The 
extra electrons occupy the 5d states, while the occupation of the 
6s and 6p states remains almost the same. 

All these charge rearrangements are caused mainly by the 
interaction of the large surface dipoles on both sides of the 
symmetrical clusters. The net result is a reduced charge transfer 
from the metal adsorbate to the substrate, and a larger fraction of 
this charge is localized near the surface. As a consequence, the 
surface dipoles are reduced significantly. The calculated surface 
dipole (in the slab-geometry case defined as ± E ± q^lz^ /2) is 
reduced for Cs adsorption from 15 to 11 Debyes, for Ba adsorption 
from 25 to 20 Debyes, and for BaO adsorption from 38 to 31 Debyes. 
For comparison, the surface dipole per BaO molecule adsorbed in 
adjacent sites on a W 14 surface cluster is reduced from 38 to 25 
Debyes. The interaction of the antiparallel surface dipoles in a 
symmetrical slab with 5 W layers is, therefore, about half as 
strong as the interaction of parallel surfaces dipoles at high 
coverage. The screening by the intermediate W layers is not very 
effective. Solid-state results for Cs or Ba(O) on W with only 5 W 
layers have therefore to be regarded with some reservations. 

TABLE X. 

Comparison of fully relativistic slab-geometry and surface -cluster 
results for Cs, Ba, and BaO on W(100). Given are populations p, 
net atmic charges q, and dipole moments p . The adsorbate -surface 
distances are: R(Cs) = 5.2, R(Ba) = 6.0, and R(O) = 1.1 a 0 . 


Property Cs/W 9 Cs 2 /W 14 Ba/W 9 Ba 2 /W 14 BaO/W 9 Ba 2 0 2 /W 14 


p (Cs/Ba, 

6sp) 

0 . 

15 

0 . 

17*2 

0. 

43 

0 . 

39*2 

0 . 

15 

0 . 

15*2 

p (Cs/Ba, 

5d) 

0 . 

24 

0 . 

33*2 

0. 

53 

0. 

72*2 

0 . 

34 

0 . 

54*2 

q (Cs/Ba) 


+ 0 . 

61 

+ 0 . 

50*2 

+1. 

04 

+ 0. 

89*2 

+ 1 . 

51 

+ 1. 

31*2 

q (0) 


, , , 

. . 

• . . 

. . 

. . . 

. . 

. . . 

. . 

-0 . 

24 

-0 . 

21*2 

q(W x ) 


+ 0 . 

04 

-0 . 

06*2 

-0. 

21 

-0. 

34*2 

+ 0 . 

06 

-0 . 

18*2 

q(W„) 


-0 . 

27 

-0 . 

30*2 

-0. 

30 

- 0 . 

29*2 

-0 . 

33 

-0 . 

32*2 

q (Wm) 


-0 . 

38 

-0 . 

14*2 

-0. 

53 

- 0 . 

26*2 

-1. 

00 

-0 . 

60*2 

p. (Debyes) 

15 

.2 

11 

.0 

25 

.1 

19 

.7 

37 

. 5 

30 

. 9 
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3.2 Work-Function/Coveraae Curves 


The important question now to be addressed is if and how the 
cluster properties for individual and interacting dipoles can be 
used to obtain information about higher coverages and possibly 
minimum work functions for actual cathode surfaces . 


3.2.1 Depolarization Model 

The decrease in the work functi on fro m 4. 5-5.5 eV for the bare 
metal substrates to around 2 eV for the cathode surfaces is caused 
by the adsorption of Ba and 0 and the creation of a surface dipole 
that lowers the electrostatic potential for electron emission. At 
small coverages, a linear variation of the work function with 
adsorbate coverage is experimentally observed. At higher 
coverages, the increasing repulsive interaction between the surface 
dipoles causes the work- function/coverage curve to deviate from the 
initially linear behavior, and a minimum is eventually reached. At 
this point the beneficial effect of an additional dipole is 
completely canceled by the destructive interaction with all other 
dipoles . 

The work function change, A $ = is given by the Helmholtz 

equation, 

A$ = - e^tn(2e 0 )' 1 , (4) 

where n is the effective surface dipole at coverage n and e 0 is the 
permittivity of vacuum. 

The effective surface dipole can be expressed as 

M = fi 0 - at , (5) 

where fi 0 is the initial surface dipole at low coverages, a an 
effective polarizability, and f the electric field resulting from 
all other dipoles. The electric field f for a simple square 
network of point-dipoles has been evaluated by Topping (1927) and 
is given by 

f = 9/m 3/2 (47re 0 ) - 1 . (6) 

Substituting f into Eqn. (5) , solving for n, and then substituting 
fi into Eqn. (4) yields 


A <f> = - 1.88 eji 0 n / (l+cotn 3/2 ) , (7) 

with /i 0 in Debyes, a in A 3 , and n in 10 1S atoms/cm 2 . The constant c 
is dependent on the actual overlayer formed and on the validity of 
the point-dipole approximation. 
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Eqn. (7) has experimentally been verified for alkali 
adsorption on single-crystal surfaces of W (Schmidt and Gomer, 
1966) and was found to be valid up to the work-function minimum. 
For larger coverages, the work function is more and more determined 
by the properties of the adsorbate itself, which are no longer 
described by the depolarization model. At sufficiently large 
coverages the work function of the pure adsorbate is eventually 
reached. 

By theoretically determining the initial dipole moments fi 0 and 
the polarizabilities a, the work-function curves should be fully 
defined up to and somewhat beyond the minimum. To test the 
validity of this theory, the two important microscopic properties, 
Mo and a, are calculated with the cluster approach. The constant 
c is determined by making an empirical adjustment for one reference 
system. 

For the initial surface dipoles ix or the calculated dipole 
moments in Table IX for the clusters with the smallest surface 
coverage are used. The polarizabilities Of are calculated from the 
induced dipole moments /x lnd (relative to fi Q ) for the clusters with 
the highest surface coverage according to 

a - Mind / Eifi , (8) 

where E^ is the electric field produced by the other dipoles in 
the cluster. The electric field f ± due to a dipole /z 0 at distance 
R is given by 

fi = Mo [ R 2 + (0 . 5[i 0 / q 0 ) 2 ] 2/3 . (9) 

Atomic units are used in evaluating Eqns. (8) and (9), i.e., 
1 a.u. (polarizability) = 0.1482 A 3 and 1 a.u. (dipole moment) = 
2.542 Debyes. R is in a 0 and q in electrons. 


3.2.2 BaO on Tungsten, Osmium, and Platinum 

The depolarization model is first applied to a comparative 
study of the work- function curves for BaO on W(100) , Os (1010) , and 
Pt(lll) . The initial surface dipoles // g and polarizabilities a are 
calculated as outlined in the previous section. The parameter c in 
Eqn. (7) has been determined such that the change in the work 
function for the model B-type Ba0/W(100) surface equals 2.65 eV, 
with 0 min =2.0 eV, which yields c = 2.7. This value is treated as 
a universal constant for simple square networks of dipoles and is 
used for all other systems with cubic crystal structures. For BaO 
on the hep crystal structure of Os, an equitriangular network is 
assumed. (The assumption of a square network is not reasonable and 
would have led to an unrealistically low work function of 1.3 eV.) 
The parameter for the equitriangular network is geometrically 
related to that for the square network by c hex = 11/9 c = 3.3, so 
that no new adjustment is needed. 
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The calculated dipole data and derived work-function changes 
A </>, minimum work functions 0 min , and optimum coverages n^ are 
presented in Table XT* Thxs table shows that the initial^ surface 
dipoles pi 0 increase in the series BaO on_ W(100), Os (1010), _and 
Pt(lll), while a shows a minimum for Os (1010). For Ba0 2 /0s (1010) , 
with Ba bridge -bonded and an excess of oxygen present in the 
cluster, both (i 0 and a are smaller than for BaO/W(100) . 


TABLE XI. 

Calculated initial surface dipoles fj. 0 (Debyes) , polarizabilities a 
(A 3 ) , initial, changes A, and minimum work functions <j> (eV) , and 
optimum coverages n„, in (in 10 15 atoms/cm 2 ) for BaO on W(100), 
Os (1010), and Pt(lll). A <f> for BaO/W(100) has been adjusted. 


System 

Mo 

a 

<t>o 

A0 

4*min 

^min 

BaO/W (100) 

30 . 0 

14.0 

4.65 

2.65 

2 . 00 

0 . 14 

BaO/Os (1010) 

40.5 

12.3 

5.2 

3.4 

1.8 

0.13 

BaO z /Os (1010) 

16.6 

9.8 

5.2 

1 . 6 

3 . 6 

0 . 16 

BaO/Pt (111) 

51.9 

27.8 

5.3 

2.9 

2.4 

0 . 09 


For the minimum work functions it is found that Ba0/0s(1010) 
has a lower and BaO/Pt(lll) a higher work function than the model 
BaO/W (100) surface. The calculated decrease in the effective work 
function for Os (0.2 eV) and the increase for Pt (0.4 eV) relative 
to W are in excellent agreement with the observed emission 
properties of Os(-W) and Pt(-W) dispenser cathodes. The work 
function for the Ba-bridge-bonded, excess-oxygen Ba0 2 /0s system is 
found to be significantly higher than for all other systems, which 
is consistent with the observed poisoning caused by excess oxygen. 
The calculated work-function curves are shown in Fig. 10. it 
should be noted that the work- function curves for coverages larger 
than n^n are spurious, because the actual work functions at higher 
coverages are known to be similar to those at the minimum or may 
even drift lower. These effects, which are due to the_ formation of 
multiple overlayers, are not included in the theoretical model. 

In conclusion, while the initial slope of the work-function 
curves is determined by fi ol the minimum work function and optimum 
coverage are dependent on the interplay between the initial dipole 
and the depolarization, and also the substrate work function (f ) 0 . 
In spite of the fact that Mo exhibits a linear behavior in the 
series BaO on W, Os, and Pt, the lowest work function and the 
largest change, A cf>, are found for BaO/Os. 
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FIGURE 10. 


Calculated work-function curves for BaO on W, Os, and Pt . 



Coverage (10 15 atoms/cm 2 ) 
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3.2.3 Other Low-Work- Function Surfaces 


As a critical test for the validity of the theoretical 
approach for the determination of work-function/coverage curves for 
low-work- function surfaces, results are now discussed for Ba and Cs 
on W(100), and BaO on W(110). The calculated dipole and work- 
function data are given in Table XXX and the work-function curves 
are shown in Fig. 11. The results for BaO/W(poly) are based on the 
assumption that /x 0 and a for the polycrystalline surface are 
similar to W(100), which appears to be justified by the fact that 
a reasonable value for the minimum work function (1.89 eV) is 
obtained . 

For Ba/W (100) , the calculated minimum work function is 2.62 
eV, which is in excellent agreement with the experimental value of 
2.63 eV (Lamouri and Krainsky, 1992). This result confirms the 
validity of the empirical adjustment for the work-function change 
of' BaO/W (100) through the choice of the parameter c in Eqn. (7) . 
The adjustment could just as well be said as having been made for 
Ba/W (100) [or BaO/W(110) , see below] . The calculated Ba binding 
energy for Ba/W(100) of 3.6 eV for the highest formal coverage 
investigated compares well the experimental value of 3 . 8 eV for Ba 
desorption from polycrystalline W (Forman, 1984) . .The calculated 
optimum coverage of 0.16 monolayers (1 monolayer = 10 15 atoms/cm ) , 
however, is significantly smaller than the expected half (0.5) 
monolayer. The small coverages appear to be a problem with the Ba- 
related systems, but the work-function values do not seem to be 
affected. 

For Cs/W (100) , the minimum work function is calculated at 
about 1.4 eV, or 0.2 eV below the experimental value of Swanson and 


TABLE XII. 

Calculated initial surface dipoles fj. Q (Debyes) , polarizabilities or 
(A 3 ), initial, changes A, and minimum work functions <f> (eV) , and 
optimum coverages (in 10 15 atoms/cm 2 ) for BaO, Ba, and Cs on 

W ( 100) , and BaO on W(110). . 


System 

Mo 

Of 

<t>o 

A<f> 

$min 

^min 

BaO/W (100) 
BaO/W (poly) 

30.0 

(30.0) 

14.0 

(14.0) 

4.65 

4.54 

2.65 

2.65 

2 . 00 
1.89 

0 . 14 
0 . 14 

Ba/W (100) 
Cs/W (100) 

20.5 

18.4 

11.8 
4 . 9 

4.65 
4 . 65 

2.03 

3.27 

2.62 

1.38 

0.16 

0.28 

BaO/W (110) 

42.2 

19.4 

5.33 

3.00 

2.33 

0 . 11 
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Strayer (1968) . In relation to the full work- function change of 
over 3 eV, this result is quite respectable in that the theory 
overshoots the total effect by only 7 per cent . The minimum in the 
work function is obtained for a Cs coverage of 0.28 monolayers, 
which is in excellent agreement with the experimental value of 0.25 
monolayers. The calculated Cs binding energy is fairly stable for 
the different clusters considered, 2. 9-3. 2 eV, and in very good 
agreement with the experimental value of 2.8 eV. 

For the BaO/W system, additional calculations have been 
carried out for different geometries of BaO on W(100) . When the 
height of Ba above the 0/W(100) surface is decreased or increased 
by 10%, the surface-dipole properties change significantly and work 
functions of 1.8 and 1.4 eV are obtained, respectively. While 
these values are not completely out of range, they show that the 
actual geometries are critical for an accurate determination of the 
work functions. 

BaO/W (110) has also been investigated. This surface is of 
particular importance because it makes up about 50% of the surface 
area of typical B-type cathodes. It is not well known, however, to 
what extent this surface (or other surfaces) contributes to the 
emission. The initial bare work function for W(llO) of 5.33 eV is 
significantly higher than that for (100) or polycrystalline W of 
4. 5-4. 6 eV. For a vertical arrangement of BaO on W(110), the 
calculated change in the work function of 3.0 eV at a coverage of 
0.1 monolayer is in excellent agreement with experiment (D. Mueller 
et al . , 1990) . The resulting work function of 2.3 eV is, however, 
too high to be relevant for cathode operation. 

The calculated Ba binding energies for the standard geometries 
of BaO/W (100) and BaO/W(110) are 5.4 and 4.6 eV, respectively, and 
compare well with the experimental value of 4.8±0.2 eV for BaO on 
polycrystalline W. 

The important issue of the contributions from different 
crystal faces to the emission of B-type cathode surfaces clearly 
requires much more work. The knowledge of accurate geometries is 
crucial, but reliable experimental information is unfortunately not 
available. Efforts are therefore planned to determine the relevant 
geometries from either a comparison of theoretical and experimental 
results or directly from first-principles calculations. With the 
good quality and consistency of the results obtained so far, the 
theoretical approach shows promise that the work functions for 
different crystal faces and alloy substrates can then be determined 
with high accuracy. 
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3.3 BaO/W Electronic Structure 


3.3.1 Densities of States 

Adsorbate interactions have been shown in the previous section 
to play an important role in determining the minimum work functions 
and optimum coverages on cathode surfaces. In order to better 
visualize and compare the results for clusters with isolated and 
interacting dipoles, total, atom-, and 1-projected densities of 
states (DOS) have been generated for BaO/W (100) from fully 
relativistic calculations. 

In Fig. 12, the densities of states (DOS) are compared for 
BaO/W 9 and Ba 2 0 2 /W 14 . The major differences are the splitting of the 
Ba 5p levels due to the Ba-Ba interaction and an increase in the 
DOS just below the cluster Fermi level for Ba 2 0 2 /W 14 . The peak 
positions of the major features are virtually unchanged. The DOS 
distributions have been generated by broadening the discrete energy 
levels of the clusters with Gaussian functions. Half -widths of 
0.30 and 0.15 eV were used for BaO/W 9 and Ba 2 0 2 /W 14 , respectively. 
The regime of tungsten-derived DOS between -5 and +5 eV is very 
well covered in both clusters, as indicated with the vertical lines 
on the horizontal axis which represent the discrete energy levels. 

In Figs. 13 and 14, atom-projected DOS are shown for one Ba, 
O, W x , W IX , and W xxx atom each in BaO/W 9 and Ba 2 0 2 /W 14 , respectively. 
In Ba 2 0 2 /W 14 , the W atoms in the center of the cluster are used. The 
comparison demonstrates that the peak positions and overall 
structure are virtually unchanged for the oxygen- and tungsten- 
related densities of states. The largest difference is observed 
for the barium- related DOS which is significantly broadened in 
Ba 2 0 2 /W 14 relative to BaO/W 9 , due to the strong Ba-Ba interaction. 

In Fig. 15, the Ba portion of the DOS and the projections on 
azimuthal quantum numbers 1=0, 1,2 are shown for both clusters on 
identical scales. The 6s and 6p (1=0 and 1) states of Ba are 
largely unoccupied, while the 5p (1=1) semi-core states are fully 
occupied. The Ba 6p states, which are somewhat higher in energy 
than the 6s states, are mixed with W states and are therefore not 
clearly seen. The 5d (1=2) states are largely unoccupied, but 
somewhat more occupied in Ba 2 0 2 /W 14 than in BaO/W 9 . The population 
analyses give Ba 5d occupations of 0.62 and 0.34 electrons per Ba 
atom in Ba 2 0 2 /W 14 and BaO/W 9 , respectively. For higher coverages the 
Ba 5d states will become more and more occupied and reduce the Ba 
oxidation level. 
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FIGURE 12. Densities of states for BaO/W 9 and Ba 2 0 2 /W 14 . 
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FIGURE 13. Projected densities of states for BaO/W 9 . 


BaO/W 9 



ENERGY RELATIVE TO THE FERMI LEVEL (eV) 


33 










36 



Contour plots are shown in Fig. 16 for the electron densities 
of BaO/W 9 and Ba 2 0 2 /W 14 in a plane perpendicular to the surface. A 
strong interaction is seen for the Ba, 0, and second-layer W atoms 
in the vertical orientation. In Ba 2 0 2 /W 14 , the barium atoms exhibit 
a large lateral interaction, while the electron density in the 0-W IX 
area is not affected. This observation correlates very well with 
the results for the respective densities of states in Figs, 13 and 
14, which show virtually no changes for the 0 and W states but a 
significant broadening of the Ba states. 


3.3.2 BaO Geometry Variation 

All previous calculations for Ba0/W(100) with Ba on top of 0 
were based on the geometry determined by Tuck and Norman et al . 
(1986 and 1987) from a fit of surface models to experimental SEXAFS 
data. Their analysis led to __R (Ba-surface) = 6.0 a 0 and R(Ba-O) = 
4.9 a 0 (or 3.2 and 2.6 A, respectively). The subsequent 
measurements of Shih et al . (1988), however, gave 2.3 A for the 
Ba-0 distance in B-type cathode surfaces. A comparison of the 
calculated binding energy for the O 2p level of 8.5 eV with the 
measured value from D. Mueller et al. (1988) for Ba0/W(100) of 5.9 
eV also hinted at a possibly smaller Ba-0 distance than as reported 
by Tuck and Norman et al. 

It is therefore of particular importance to investigate the 
effects of different geometries on the charge transfer, the dipole 
moment, and the 0 2p binding energy in the BaO/W surface complex. 
Calculations have been carried out for BaO/W 9 assuming the vertical 
adsorbate orientation, with Ba-surface distances of 6.0, 6.3, and 
6.6 a 0 and the two proposed Ba-0 distances of 4.9 and 4.3 a 0 (2.6 
and 2.3 A) . 


TABLE XIII. 

Calculated net atomic charges q (electrons) , dipole moments fj. 
(Debyes) , and Ba dipole energies E Ba (eV) for different geometries 
of BaO/W 9 at the fully relativistic level. R (Ba) = R (Ba-surface) 
and R(Ba-O) are given in Bohrs (1 a 0 = 0.5292 A) . 


R(Ba,Ba-0) 

^Ba 

q 0 


^Ba 

6.0, 4.9 

1.51 

-0.24 

37.5 

6.4 

6.3, 4.9 

1.51 

-0.32 

38.0 

6.3 

6.6, 4.9 

1.49 

-0.39 

38.2 

6 . 0 

6.0, 4.3 

1.36 

-0.36 

33.4 

5.2 

6.3, 4.3 

1.33 

-0.43 

32.6 

5 . 0 

6.6, 4.3 

1.29 

-0.48 

31.9 

4 . 7 
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The results for the surface -dipole properties and the 
electronic binding energies are presented in Tables XIII and XIV, 
respectively. 

The following effects are observed: When R(Ba) = R (Ba-surface) 
is increased or R(Ba-O) is decreased, the net charge on 0 increases 
which leads to smaller 0 2p binding energies. When R(Ba-O) is 
changed from 4.9 to 4.3 a 0 , the net charge on Ba, the cluster 
dipole moment, and the Ba desorption energy all decrease. The 
decrease in the net positive charge on Ba leads to slightly smaller 
electronic Ba binding energies, while the W electronic levels are 
not affected. The increase in the negative charge on 0 for larger 
Ba-surface distances gives rise to increases in the negative 0-W 
contributions to the dipole moment, which offset the increases due 
to the larger Ba-surface distances and leave the total dipole 
moments within R{Ba-0) = 4.9 and 4,3 a 0 virtually unchanged. 


A comparison of the results for the original geometry with 
R(Ba,Ba-0) = 6.0, 4.9 a 0 to those based on 6.6, 4.3 a 0 shows 

significant improvements with respect to the experimental data. 


TABLE XIV. 

Fully relativistically calculated electronic binding energies (in 
eV) relative to the highest occupied cluster level for different 
geometries of BaO/W 9 . The final results for R(Ba,Ba-0) = 6.6, 4.3 
a„ include the final-state contributions (AFS) . R(Ba) 

R (Ba-surface) and R(Ba-O) are given in Bohrs (1 a 0 = 0.5292 A) . 


R (Ba, Ba-O) 

^Pl/2 

W, 

^P3/2 

4 f 5/2 

7/2 

5s 

-Ba-- 

5Pl/2 

5p3/2 

2s 

0 

2pcr 

2pir 

w„ 

6s 

6.0, 4.9 

45.8 

36.0 

30.0 

27.7 

30.9 

17.3 

15.2 

22 . 7 

8.6 

8 . 5 

5 . 6 

6.3, 4.9 

45.8 

36.0 

30.0 

27.8 

30.8 

17.2 

15.1 

21.8 

7.8 

7.7 

5 . 6 

6.6, 4.9 

45.8 

36.0 

30.1 

27.8 

30.6 

17.0 

14.9 

21 . 0 

7.3 

6 . 9 

5.6 

6.0, 4.3 

45 . 8 

36.0 

30.0 

27.8 

29.2 

15.6 

13.6 

21.6 

7.3 

7.0 

5 . 6 

6.3, 4.3 

45.8 

36.0 

30.1 

27.8 

28 . 9 

15.3 

13.2 

20.8 

6.9 

6.2 

5 . 5 

6.6, 4.3 

45.8 

36.0 

30.1 

27.8 

28.6 

15.0 

12 . 9 

20.1 

6.6 

5 . 5 

5 . 3 

AFS 

0 . 9 

0.7 

4.1 

4 . 0 

3.1 

2.3 

2.5 

1.4 

0 . 1 

0.7 

0.2 

Final 

46.7 

36.7 

34.2 

31.8 

31.7 

17.3 

15.4 

21.5 

6.7 

6.2 

5 . 5 

Expt . a 

.... 

36.8 

33.5 

31.3 

.... 

16.6 

14.6 

21.9 


5.9 


AExpt . 

.... 

-0.1 

+ 0.7 

+ 0.5 

.... 

+ 0.7 

+0.8 

-0.4 

+ 0 . 

, 8/-G 

1.4 


a D. Mueller et al . (1988) 
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The calculated Ba dipole (desorption) energy decreases from 6.4 to 
4.7 eV (expt. 4.8 eV) and the 0 2p binding energy decreases from 
~8.5 to around 6 eV (expt. 5.9 eV) . Because the W IX 6s state has 
a similar binding energy, a substantial mixing of the W XI 6s and 0 
2pa orbitals occurs, which increases with decreasing energy 
separation between these two orbitals. For R(Ba,Ba-0) = 6.6, 4.3 
a 0 , the state at 6.6 eV actually contains a larger contribution 
from W xx 6s (25%) than from 0 2pa (15%), and the 0 2pa character is 
found mainly in the state at 5.3 eV. 

For the final geometry, R (Ba, Ba-O) = 6.6, 4.3 a 0 , additional 
calculations have been carried out in order to investigate 
corrections to the orbital energies of the N-electron system due to 
the reorganization of the charge distribution upon ionization. 
This is done by performing so-called transition-state calculations 
in which 1/2 electron is removed from the state under 
consideration. The (N-l/2) -electron system is the "transition 
state" between the N- and (N-l) -electron systems. This procedure 
is known to provide good estimates for the differential 
reorganization, or final-state, effects (AFS) . A separate 
calculation is required for each state, and the electron binding 
energies are now measured relative to the transition-state energy 
of the highest occupied orbital . 

In comparison with photoemission spectra, only the differences 
in the calculated energy shifts are relevant in cluster 
calculations. The differential effect for the electronic state at 
the cluster Fermi level is zero, by definition, since it is used as 
the reference. The results in Table XIV show that the final-state 
contributions (AFS) increase with increasing binding energy and 
localization. For the W xx 6s and 0 2pa states, which are 
significantly mixed, the correction is only 0.1 to 0.2 eV. AFS is 
of the order of 1 eV for the more localized O 2p7r, 0 2s, and W 5p 
electronic states, 2-3 eV for Ba 5p and 5s, and 4 eV for W 4f. 

After inclusion of these final-state corrections, all 
electronic binding energies below 50 eV for BaO/W 9 with R (Ba, Ba-O) 
= 6.6, 4.3 a 0 are within 1 eV of the photoemission energies from D. 
Mueller et al . (1988), while the results for the Ba- and O-related 
energies with R(Ba,Ba-0) = 6.0, 4.9 a 0 would have been too high by 
2 eV or more. Fig. 17 shows the densities of states for these two 
geometries and the final results for R(Ba,Ba-0) = 6.6, 4.3 a 0 . The 
corresponding experimental peak positions are indicated by arrows. 
Interestingly, the Ba binding energy for the BaO/W 9 cluster with 
the modified geometry, 4.7 eV, is within the margin of experimental 
uncertainty (4.8±0.2 eV) . 

The good agreement between theory and experiment for the 
modified geometry is strong evidence that the distances proposed by 
Tuck and Norman et al . for Ba0/W(100) are most likely not correct. 
The present calculations clearly favor the smaller Ba-O distance 
determined by Shih et al . (1988), although in his model Ba is 
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bridge-bonded to 0, and no fine structure indicative of a Ba-W 
distance was observed. The better agreement for a larger 
Ba-surface distance than as proposed by Tuck et al. (1986) could 
also be a result of the absence of adsorbate interactions in 
BaO/Wg, which, from the results for Ba 2 0 2 /W 14 , can be expected to 
lead to a somewhat smaller charge transfer, smaller binding 
energies, and possibly similarly good agreement with experimental 
binding energies at a slightly smaller Ba-surface distance. 


FIGURE 17 . 

Density of states for Ba0/W 9 with R(Ba, Ba-0) = 6.0, 4.9 a 0 
(dotted) and with R(Ba, Ba-0) = 6.6, 4.3 a 0 (dashed, 
and full line after inclusion of final-state effects) . 

The arrows indicate experimental peak positions. 
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Ill .litJkiii 


The general conclusions with respect to the charge transfer 
and surface dipole, on the other hand, are qualitatively similar 
for the modified and original geometries: Ba is strongly oxidized 
with q Ba = 1.3 electrons (1.5 orig.) and 0 is only slightly reduced 
with q 0 = -0.5 electrons (-0.24 orig.), which gives rise to a very 
large surface dipole of 32 Debyes (37 orig.). 

For the determination of the final-state effects for the 5p 
and 4f core states of the four equivalent W surface atoms in BaO/W 9 
the following approximation has been employed: A total of 1/2 
electron was removed from the symmetry- equivalent state, which 
means that only 1/8 electron was removed per site. The calculated 
energy shift from the transition-state calculation is then scaled 
by a factor of four. 

Because the W r atoms are related by symmetry, a rigorous 
transition-state calculation with 1/2 electron removed from one W x 
core state would reduce the fourfold symmetry, C 4v/ to C s with just 
one reflection plane remaining. This would substantially increase 
the computational effort, in particular, since a separate 
calculation is required for each single state. 

The validity of the above approximation is demonstrated by a 
comparison of the results for the W x and W IX core states in Table 
XV. The final -state contribution for the W x 5p and 4f states is 
calculated to be 0.2 and 1.0 eV, respectively, or after scaling 
close to 1 and 4 eV. The latter values compare very well with the 
directly calculated contributions for W„ where no scaling is 
needed . 


TABLE XV. 

Orbital energies Ae, relative to the cluster Fermi level, final - 
state contributions AFS, final calculated, and experimental binding 
energies for tungsten-related states of BaO/W(100), in eV. 
R (Ba-surface) = 6.6 a 0 and R(Ba-O) = 4.3 a 0 ~2.3 A on W 9 . 


s— * 

State 

-Ae 

AFS AFS (scaled) Final 

Expt . 

AExpt . 

H 

5 Pi/2 

45 . 8 

0.2 

0 . 9 

46.7 






5p 3 /2 

36 . 0 

0.2 

0.7 

36.7 

36 . 8 

-0 . 1 


4 f 3/2 

30.1 

1 . 0 

4 . 1 

34.2 

33.5 

+ 0 . 7 

t 4 

4f 5 / 2 

27.8 

1.0 

4 . 0 

31.8 

31.3 

+ 0.5 

F-'J 

*P— 

^ii 5p x / 2 

46 . 0 

1 . 1 

1 . 1 

47.1 

— 

— 

- . 

5p 3 /2 

36.5 

1.0 

1 . 0 

37.5 

36 . 8 

+ 0.7 

- 7 

4 f 3/2 

29.8 

4.5 

4.5 

34.3 

33.5 

+ 0.8 


4 f 5/2 

27.6 

4.3 

4.3 

31.9 

31.3 

+ 0 . 6 
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3.3.3 Effects of Bulk and Surface Environments 


a first stGp toward the simulation of bulk and surface 
environments in atomic cluster calculations, W 9 and BaO/W 9 clusters 
have been investigated in which the sphere radii for the various 
atoms have been modified to more realistically represent the 
conditions that are present in real, extended systems. 

The standard atomic sphere radii for an isolated cluster are 
chosen as 88 percent of the radii determined from the superposition 
of nonrelativist ic atomic charge densities. A standard ^ body- 
centered-cubic W 9 cluster, for example, has eight identical W 
surface atoms and one inner W atom. 

If used as a surface model, however, the third-layer tungsten 
atoms, W m , should be treated differently from W If since they have 
tungsten neighbors in an extended system with similar 
nearest-neighbor distances as W„. A set of modified calculations 
has therefore been carried out for W 9 and BaO/W 9 , in which the 
sphere radii for the four W IXI atoms were reduced to those of W„. 
In a third set of calculations, the sphere radii were taken from 
more truly embedded systems, W 34 and Ba 5 0 5 /W 34 , which leads to a 
reduction of all sphere radii relative to the standard clusters . 
W 34 and Ba 5 0 5 /W 34 are derived from W 9 and BaO/W 9 by adding the eight 
nearest neighbors in the first and third W layers, the four nearest 
neighbors in the second layer, a fourth W layer similar to the 
second layer, and in the adsorbate cluster the four nearest BaO 
molecules . 

Fully relativistic results for the series of calculations for 
W 9 are given in Table XVI and Figs. 18-20. The charge distribution 
is symmetrical for standard W 9 , with about 0.05 electrons 
transferred from each of the surface atoms to the central W xx atom. 
For W 9 with modified and embedded radii, the W XI atom acquires 0.7 
to 1.0 additional electrons, mainly from the W IIX atoms, . This 
charge redistribution sets up negative dipole moments perpendicular 
to the Wj "surface" layer with excess negative charge at or near 
the surface as present on real metal surfaces. The reduction of 
the sphere radii in the series leads to a stabilization of the 
electronic cluster levels, as expressed by the shift in the highest 
occupied state, £max' which is indicative of an increase in bulk 
character. 

In Figs. 18, 19, and 20 the density of states for W 9 is 
decomposed into contributions arising from the surface, W(S) , 
subsurface, W(S-l), and third- layer, W(S-2), atoms, respectively. 
A comparison with the projected densities of states (see insets) 
from the solid-state calculation of Wimmer et al.(1983) for a 
5_layer W film shows good agreement for the surface and subsurface 
layers . 
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TABLE XVI. 


Net charges q (in electrons) , dipole moments /z, and energies e max of 
the highest occupied state for W 9 clusters with standard, modif Ted, 
and atomic sphere radii taken from an embedded cluster (W 34 ) 


Property 

stand, radii 3 

mod. radii b 

emb. radii c 

q(W x ) 

+ 0.19 

-0.44 

+0.30 

q(w IX ) 

-0.38 

-1.08 

-1.39 

q(w XII ) 

+ 0.19 

+1.52 

+ 1.09 

/z (Debyes) 

0.00 

-8.16 

-4.06 

Emax (eV) 

-5.07 

-5.78 

-6.15 


* R(W I ) =3 . 161, R (W xx ) =2 . 672 , R (W xxx ) =3 . 161, R(outer sphere) =8 . 340 a 0 . 
15 As for a, except R(w m ) = R(W XX ). 

c R(W X )=2.850, R (W IX ) =2 . 593 , R (W xxx ) =2 . 782 , R(outer sphere) =8 . 340 a 0 . 


The densities of states in Fig. 18 for the surface atoms (S) 
of W, with different atomic sphere radii show the prominent surface 
state at -1 eV, four peaks between -2 and -6 eV, the broad density 
of unoccupied states around +1 eV, and additional features at 
higher energies. The surface structure for this simple cluster is 
in remarkably good agreement with the (non-relativistic) 
solid-state result for the 5-layer W(100) film which is shown on 
the inset in the middle figure. The intensity distribution is 
reasonably well reproduced by all clusters, while the width of the 
occupied DOS is better described by W 9 with the radii taken from 
the embedded cluster. 

The subsurface-atom (S-l) densities of states for the 
different W 9 clusters are compared in Fig. 19. The surface state 
and the high DOS around the Fermi level are now completely missing. 
Higher DOS are found instead below -2 and above +1 eV, with the 
consequence that the bulk- like DOS exhibits a minimum where the 
surface DOS has its maximum. The sharply peaked DOS is character- 
istic of d-electron metals. Again, the simple cluster result is in 
good general agreement with the solid-state calculation. 

The DOS for the third-layer atoms (S-2) in Fig. 20 is for 
standard W 9 identical and for the other two calculations similar to 
the respective surface DOS of Fig. 18. Despite the large changes 
observed in the calculated charge distributions, the projected 
densities of states are fairly similar for the different members of 
the series, and only the second- layer W can be considered as having 
bulk-like properties. 
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FIGURE 18. 


Densities of states for the surface atoms (S) of W 9 with 
standard, modified, and atomic sphere radii taken from 

an embedded cluster. 
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FIGURE 19. 

Densities of states for the subsurface atoms (S-l) of W 9 with 
standard, modified, and atomic sphere radii taken from an 

embedded cluster. 
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FIGURE 20. 
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The reduction of the sphere radii in the BaO/W 9 cluster with 
R(Ba,Ba-0) = 6.0, 4.9 a 0 has the following effects on the fully 
relativistically calculated charge distribution and dipole 
properties, as shown in Table XVII: The charge on Ba increases 
steadily in the series from 1.51 to 1.62 and 1.83 electrons, which 
is likely to be related to the stabilization of the tungsten 
electronic structure already observed for bare W 9 . The screened 
negative image charge is now found mainly in the second rather than 
in the third W layer. A comparison with the charge distributions 
in the respective bare W 9 clusters (Table XVI) shows, ‘however , that 
the Wjj atom already had a similarly large negative charge such that 
the actual charge transfer is still into the third tungsten layer. 
The dipole moments, therefore, remain similarly large, while the Ba 
binding energy increases substantially due to the increasing 
positive charge on Ba. It is also worth noting that the reduction 
of the W x sphere radius from 3.16 to 2.85 a 0 in BaO/W 9 with embedded 
radii leads to increases in the oxidation level of the W surface 
layer and the negative charge on oxygen. 


TABLE XVII. 

Net charges q (in electrons) , dipole moments /i, Ba dipole energies 
E^, and energies e ^ of the highest occupied state for BaO/W 9 
clusters with standard, modified, and atomic sphere radii taken 
from an embedded cluster (Ba 5 0 5 /W 34 ) . R(Ba-O) = 4.9 a 0 ~2 . 6 A. 


Property 

stand. radii a 

mod. radii b 

emb. radii 0 

q (Ba) 

+1.51 

+1.62 

+ 1.83 

q (O) 

-0.24 

-0.23 

-0.73 

q (Substrate) 

-1.27 

-1.39 

-1 . 10 

q (W r ) 

+0.06 

-0.30 

+ 0.73 

q(W IX ) 

-0.33 

-1.13 

-1.34 

q(w Ixl ) 

-1.00 

+0.04 

-0.49 

fi (Debyes) 

37.5 

37.2 

43.2 

Eea (eV) 

6.4 

7.9 

9.8 

(eV) 

-5.39 

-5.95 

-6.79 


a R(Ba) =3 . 769, R(0)=1.777, R(W I )=3.000, R (W ri ) =2 . 613 , R (W TII ) =3 . 160 , 
R (outer sphere) =10 . 4 86 a 0 . 
b As for a, except R(W IIX ) = R(W IX ) . 

c R (Ba) =3.323, R(0)=1.763, R(W X )=2.684, R(W„) =2 . 544 , R (W in ) =2 . 673 , 
R (outer sphere) =10 . 486 a 0 . 
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The densities of states for BaO/W 9 with standard, modified, 
and atomic sphere radii from the embedded cluster are compared in 
Fig. 21. The chief differences are (i) changes in the W-related 
DOS outside the regime of high surface DOS around the cluster Fermi 
level and (ii) a shift of the Ba 5p states at -15 eV by 1 . 5 eV to 
smaller binding energies relative to standard BaO/W 9 . 
Surprisingly, the Ba 5p levels are not affected by the large 
reduction of the Ba sphere radius when comparing BaO/W 9 with 
modified and embedded radii. The 0 2s and 2p binding energies 
remain unchanged for all clusters. 

The results for the BaO/W electronic structure can be 
summarized as follows: 

The 0 2s and 2p binding energies can be brought in agreement 
with the photoemission data of D. Mueller et al . (1988) by 
changing the Ba-0 distance in the vertical arrangement from 
the value of 2.6 A proposed by Tuck and Norman et al . (1986 
and 1987) to the value of 2.3 A proposed by Shih et al . 
(1988) . Final-state effects are important for the 0 2s (1.4 
eV) and 0 2p tt (0.7 eV) binding energies. Adsorbate 
interactions and simulated environments do not affect the 
oxygen electronic levels. 

The Ba 5p energies also improve with a shorter Ba-0 distance. 
The final-state effects amount to 3.1 eV for Ba 5s and -2.4 eV 
for 5p. These states are strongly affected when the smaller 
sphere radii of embedded W atoms are used, but the effect of 
adsorbate interactions is small. 

All calculations so far have confirmed a large charge transfer 
from Ba to the substrate, a relatively small negative charge 
on 0, and a resulting very large surface dipole. The 
calculated Ba binding energies are in good agreement with 
experiment . 
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FIGURE 21. 


Densities of states for BaO/W 9 with standard, modified, and 
atomic sphere radii taken from an embedded cluster (Ba 5 0 5 /W 34 ) 
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4. Q/W ELECTRONIC STRUCTURE 


For the 0/W(100) system, e xten ded cluster calculations have 
been carried out and a comparison has been made with oxygen- induced 
features in angularly- resolved inverse photoemission spectra 
(Krainsky, 1988) . 

Figure 22 shows the fully relativistically calculated 
densities of states (DOS) for a 0 2 /W 14 cluster, the W 14 substrate, 
the oxygen -related DOS in 0 2 /W 14 , the difference 0 2 /W 14 -W 14 , and the 
experimentally observed spectrum of oxygen- induced features for an 
angle of incidence of 52°, which is representative of the spectra 
taken between 40° and 76°. The density of the oxygen- induced states 
in this range is in remarkably good agreement with experiment . The 
intensities and energies of the observed peaks are well reproduced 
and significantly improved from earlier results for single- 
adsorbate O/W clusters. The oxygen -related DOS exhibits only one 
major peak, the O 2p peak, at 6.3 eV below the cluster Fermi level. 
This shows that all other features in the difference spectrum 
originate from changes in the electronic structure of W due to the 
adsorption of O, and not from unoccupied O levels themselves. 

For the purpose of comparison with the above results, 
quasirelativistic calculations have been carried out for the same 
system. This comparison revealed significant differences. At the 
quasirelativistic level, the occupied O 2p states are shifted from 
-6.3 eV to -4.0 eV, and the peaks near 1 and 2 eV above E F are now 
larger in intensity than the state near 4 eV. 

Preliminary surveys have also been carried out for oxygen 
adsorped in a bridge position at room temperature and for a 
reconstructed surface at elevated temperatures by using the 
quasirelativistic approach. The adsorption of oxygen in a bridge 
position on unreconstructed W(100) is modeled with a 0 4 /W 12 cluster, 
where the substrate cluster contains two layers of six W atoms 
each. A 0 2 /W 12 cluster is used for the high- temperature case, where 
the middle row of the first -layer W atoms has been moved to the 
third layer and oxygen is now adsorbed on top of the exposed 
second-layer W atoms. For the unreconstructed surface, the results 
for 0 4 /W 12 -W 12 with oxygen at a bridge site are similar to those with 
oxygen at the fourfold-hollow site. For the reconstructed surface, 
the DOS for 0 2 /W 12 -W 12 shows a further increase in the intensity for 
the peak near 1 eV, so that none of these results provides an 
improvement in the description of the observed spectrum. 

Because of the shortcomings of the quasirelativistic results 
with respect to the intensities of the oxygen- induced unoccupied 
states, it has been decided to further explore the O/W system with 
fully relativistic calculations for larger clusters, in which 
proper surface equivalences will taken into account . 
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Difference Density of states (arb. units) 


FIGURE 22. 

Theoretical and experimental oxygen -induced features 
in the density of states for 0/W(100) . 
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5. AB INITIO BaO RESULTS 


Ab initio configuration interaction (Cl) calculations, which 
are based on self -consistent-field (SCF) or multi-configuration 
(MC) wave functions, have been carried out for BaO. The electronic 
structure on Ba is described by using a previously developed 
2 -electron pseudopotential which includes relativistic and core 
polarization (CPP) effects (Muller, 1988 ) . In Table XVIII the Cl 
and MC-CI results are compared with earlier SCF and MC-SCF data. 

The spectroscopic constants in this table include the 
equilibrium separation R e , the rotational constant B e , the harmonic 
vibrational frequency co e , the dissociation energies D e and D 0 , the 
dipole moment pt e , and the net charge on Ba q^. The calculation 
denoted as MC-CI+CPP' has been corrected for the Ba 2+ (core) - O' 
(2p a" 1 ) exchange repulsion as obtained from the superposition of the 
O' Cl orbitals with an all-electron Ba 2+ core, while for all other 
cases the O' SCF orbitals have been used. 


TABLE XVIII. 

Spectroscopic constants for BaO at the SCF, MC and Cl levels of 
calculation and from experiment (Huber and Herzberg, 1979) . 


Method 

Re ( A ) 

B e (cm* 1 ) 

co e (cm' 1 ) 

D e (eV) 

D 0 (eV) 

Me ( D ) 

qe a (e) 

SCF 

1.90 

0.325 

972 

5.9 

3 . 7 

- 12.2 

1.2 

SCF+CPP 

1.89 

0.328 

976 

6.5 

4.3 

- 11.0 

1.1 

MC-SCF+CPP 

1.92 

0.319 

813 

7.0 

4.8 

- 9.8 

1 . 0 

CI+CPP 

1.91 

0.324 

944 

8.5 

6.3 

- 10.1 

1.0 

MC-CI+CPP 

1.88 

0.331 

950 

8.9 

6.8 

- 9.2 

0 . 9 

MC-CI+CPP' 

1.91 

0.324 

944 

8.6 

6.5 

- 9.2 

0 . 9 

Expt . 

1.94 

0.313 

670 

7.8 

5.8 

- 8.0 

0 . 9 

The general 

conclusions are 

as follows: 




R e and B e are very well 

described at all 

levels 

* 



o> e is still significantly too large at the Cl level. This is 
most likely related to the approximate treatment of the Ba 2+ 
core - O' repulsion, in which the transfer of 1 electron is 
assumed for all R around R e . A more rigorous treatment would 
require the explicit consideration of the next inner shell on 
Ba . 
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The calculated dissociation energies are somewhat too small at 
the SCF and MC-SCF levels and appear to be too large at the Cl 
level. A rigorous experimental value, however, is not known. 

The calculated dipole moment and charge transfer improve 
systematically with increasing level of computational effort. 


In summary, it is found that the geometry and bonding in BaO 
is reasonably well described with the Ba 2-electron pseudopotential 
at the intermediate MC-SCF+CPP level, while a better description of 
the shape of the potential curve requires the explicit 
consideration of the next inner electronic shell on Ba and/or 
possibly also a larger basis set on oxygen. 
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The electronic structure of the high-temperature super- 
conductor YBa 2 Cu 3 0 7 has been investigated with the relativistic 
cluster approach (Muller, Physica C, 1989) . The results described 
here represent the first calculations for this system which include 
(a) all relativistic electronic effects and (b) final-state energy 
corrections for comparison with photoemission spectra. 

The atomic clusters used for YBa 2 Cu 3 0 7 are shown in Figure 23. 
These include two small clusters modeling the Cu-0 chains and 
planes, and larger clusters for two different unit cells with 3 and 
4 Cu-0 layers. 


FIGURE 23. Cluster models for YBa 2 Cu 3 0 7 . 


M tfrf 

BaCiuOe YBaOuOe 
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The valence densities of states (DOS) for the two small 
clusters and the Y 2 Ba 2 Cu 12 0 18 unit -cell cluster are compared in Fig. 
24 . The DOS is generated from the calculated energy levels which 
are indicated as small vertical lines on the horizontal axis. The 
peaked and smooth density distributions were obtained by broadening 
the cluster energy levels with Gaussians using width parameters of 
0.05 and 0.50 eV, respectively. The BaCu 4 0 6 cluster, modeling the 
chains, already shows the characteristic double peak below the 
Fermi level that is observed in photoemission spectra. For 
YBaCu 4 O e , which models the planes, additional states appear in the 
center of the double-peak structure. For the unit-cell cluster, 
the width and shape of the double peak is in good agreement with 
experiment. The additional feature at -5 eV results from the 
planes. Above the Fermi level, the states at lower energy are due 
to Y and the states at higher energy are due to Ba. 

The calculated total and partial densities of states for the 
Cu-0 valence band of Y 2 Ba 2 Cu 12 0 18 are shown in Figure 25, which shows 
the typical double maximum and additional features at higher and 
lower binding energies . The partial DOS from the cluster shows 
major copper contributions near the center of the double maximum 
and oxygen contributions more to its sides. These effects are 
analyzed in more detail in Fig. 26. In this figure, photoemission 
spectra for different photon energies are compared with calculated 
DOS for Cu-0 in which the Cu contribution is varied from 0 to 100%. 
It can be seen that the double-peak structure of the high-energy, 
hv = 70 eV, spectrum is well reproduced by the calculated Cu-0 
(both 100%) DOS, while the three -peak structure of the low-energy, 
hv = 20 eV, spectrum resembles the Cu (0%) -O (100%) DOS. Low photon 
energies, therefore, emphasize the oxygen contribution to the Cu-0 
valence band. 

In Fig. 27, the DOS are compared for different unit-cell 
clusters. [Y 2 Ba 2 Cu 12 0 18 2 "] 2+ contains two additional electrons and is 
surrounded by a "Watson" sphere with charge 2+, which represents an 
attempt to account for effects of the ionic environment in the 
crystalline material. If formal charges of +3 and +2 are assigned 
to Y and Ba, respectively, and -2 to O, then the formal charge of 
Cu in YBa 2 Cu 3 0 7 is +2.33. Because the clusters do not have the 
correct stoichiometry, the cluster of atoms in, e.g., Y 2 Ba 2 Cu 12 0 18 
requires 2 additional electrons for neutrality of the formal 
charges, which are balanced with 2 positive charges in the Watson 
sphere. Fig. 27 shows that this modification leads to a slight 
shift in the DOS of the order of -0.5 eV because of the occupation 
of a previously unoccupied orbital in Y 2 Ba 2 Cu 12 0 18 . The choice of a 
different unit cell, YBa 2 Cu 16 O 20 , with the Cu-0 planes and Y at the 
center of the cluster (instead of the Cu-0 chains and Ba) , has a 
somewhat larger effect on the DOS, but the calculated general 
structure remains the same . 







FIGURE 25. 


Total (per YBa 2 Cu 3 0 7 cell) and partial (atomic) densities 
of states as derived from a Y 2 Ba 2 Cu 12 0 18 cluster model. 
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Intensity (arb, units) 


FIGURE 26. 


Analyses of experimental and theoretical 
densities of states for Y-Ba-Cu-0 near the Fermi level. 


EXPERIMENT THEORY 




58 








Overall, the cluster DOS is in good agreement with solid-state 
calculations and photoemission data. Relative to the results from 
both calculations, the observed valence band is shifted, however, 
by about 1.5 eV to higher binding energies (Johnson et al . , 1987) . 
This shift, as discussed below, has to be attributed to true 
electron correlation effects which are not included in calculations 
that are based on the density- functional approach. 

The calculated charge distributions for the various clusters 
are given in Table XIX. For Ba, Cu, and 0 net charges of +/- 0.6 
to 0.8 electrons are obtained. The net charge for Y shows the 
largest variation, with between +0.7 and +1.3 electrons. Although 
the real charges in the bulk material are certainly not as large as 
the formal charges given above, the true net charges in the bulk 
material are likely to be larger than as obtained in the cluster 
calculations because the long-range Madelung contribution to the 
potential is missing. This contribution is expected to stabilize 
additional electronic charge on oxygen and to increase the ionicity 
of the material. Because, however, the ionic environment is 
missing near the surface of a superconductor, the cluster results 
are believed to be representative of the surface electronic 
structure . 


TABLE XIX. 


Net atomic charges (in electrons) for different Y-Ba-Cu-0 
superconductor clusters from fully relativistic calculations. 


Cluster 

Y 

Ba 

Cul 

Cu2 

01-4 

BaCu 4 O e 

• • • • 

+0.8 

+0.8 

.... 

-0.7 

YBaCu 4 0 8 

+ 1.3 

+0.8 

.... 

+0.8 

-0.6/-0.7 

Y 2 Ba 2 Cu 12 0 18 

+1.1 

+ 0.6 

+ 0.6 

+0.8 

-0 . 6/-0 . 8 

[Y 2 Ba 2 Cu 12 0 18 2 ] 

+ 1.0 

+ 0.5 

+0.6 

+0.7 

-0.7/-0.8 

YBa 2 Cu 16 O 20 

+ 0.7 

+ 0.6 

+0.7 

+ 0.8 

-0 . 6/-0 . 8 


Details of the calculated valence charge distributions for 
Y 2 Ba 2 Cu 12 0 18 and YBa 2 Cu u O 20 are compared in Table XX. Although there 
are differences in the net atomic charges, in particular for Y 
which is either at the upper and lower ends or in the center of the 
cluster (see Fig. 23) , the s, p, and d orbital populations for the 
various atoms are fairly similar in the two cases. 
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Barium has a 6sp population of 0.4 electrons and a 5d 
population of about 1 electron, i.e., a net positive charge of 
about 0.6 electrons, and oxygen has an excess of 0.6 to 0.8 
electrons in the 2p shell. The Cu 4sp population varies between 
0.6 and 1.1 electrons and the Cu 3d population between 9.3 and 9.8 
electrons. The net charge on Cu is on the order of +0.7 electrons, 
and for Cu(l) approximately 0.15 electrons smaller than for Cu(2) . 
The results for both clusters show that holes are present mainly in 
the O 2p shell. These holes are believed to be the charge carriers 
in the high- temperature superconductors. The charge distribution 
also indicates that there exist significant covalent contributions 
in the interatomic bonding near the surface. 

The Cu 3d occupations of 9. 3 -9. 6 electrons are in very good 
agreement with the estimate of 9.45 from the Larsson-Sawatzky model 
(Fujimori et al., 1987). The calculated Cu 3d spin-orbit 
splittings are 0.3 eV. 


TABLE XX. 

Valence charge distribution for unit-cell clusters of YBa 2 Cu 3 0 7 from 
relativistic Xa calculations. 


Cluster 

s 

P 

d 

sum 

net 

Y 2 Ba 2 Cu 12 0 18 





r • 

Y 

0.26 

0.24 

1.41 

1.91 

+ 1.09 

Ba 

0.17 

0.27 

1.00 

1.44 

+ 0 . 56 

Cul 

0.27 

0.81 

9.28 

10.36 

+ 0 . 64 

Cu2 

0.30 

0.30 

9.61 

10.21 

+0.79 

Ol 

2.03 

4.72 

.... 

6.75 

-0 . 75 

02 

2 . 03 

4.68 

, , , , 

6.71 

-0 . 71 

03 

2.03 

4.67 

, , , , 

6.70 

-0.70 

04 

2.06 

4.57 

.... 

6.63 

-0.63 

YBa 2 Cu 16 O 20 






Y 

0.27 

0.41 

1.60 

2.28 

+ 0.72 

Ba 

0 . 17 

0.28 

0.95 

1.40 

+ 0 . 60 

Cul 

0.27 

0.55 

9.50 

10.32 

+ 0 . 68 

Cu2 

0.31 

0.28 

9.60 

10.19 

+ 0 . 81 

01 

2 . 02 

4.68 

.... 

6.70 

-0.70 

02 

2 . 04 

4.74 

, , ( , 

6.78 

-0.78 

03 

2 . 03 

4.74 

.... 

6.77 

-0 . 77 

04 

2 . 08 

4.52 

.... 

6.60 

-0.60 


61 





Energy shifts due to the relaxation of the charge distribution 
upon ionization, which are the final-state contributions in the 
photoemission process, h ave be en investigat ed f or the first time. 
It is found that these effects are small for the valence states 
(for localized Cu 3d states up to 0.2 eV) , and that the valence- 
band width is not affected. True correlation effects have, 
therefore, to be made responsible for the observed shifts. 

The final -state contributions are, however, very large for the 
core states. They range from 160 eV for Cu Is a nd 30 eV for O Is 
to 1 eV for Cu 3s, 3p, and O 2s. The resulting core-electron 
binding energies for Y 2 Ba 2 Cu 12 0 18 are given in Table XXI . The 
calculated values are in very good agreement with experimental 
data, with deviations of 1-5 eV only. Without consideration of the 
final -state contributions, the errors would be up to one order of 
magnitude larger. 


TABLE XXI. 


Calculated final-state contributions (AFS) and core-electron 
binding energies (eV) for Y 2 Ba 2 Cu 12 0 18 . 


Atom 

State 

AFS 

Calc. 

Expt . 

Cu(l, 2) 

Is 

160 

8999, 

9004 




2s 

17 

1077, 

1077 

— 


3s 

1 

113, 

112 

— 


2Pl/2 

23 

953, 

951 

953 


3Pl/2 

1 

73, 

72 

76 


2p3/2 

23 

932, 

931 

933 


3P3/2 

1 

70, 

69 

74 

0(1,4) 

is 

30 

533, 

533 

529 


2s 

1 

18, 

17 







(531,533) 

0(2,3) 

Is 

30 

534, 

534 



2s 

1 

19 , 

19 

20 

Ba 

3d 3/2 

23 

794 

— 


3d s ; 2 

23 

778 

779 

Y 

3d 3 / 2 

10 

160 

158 


3d 5/2 

10 

158 

156 


A close examination of the results reveals that the calculated 
Cu levels are slightly lower and the 0 levels slightly higher than 
the corresponding experimental values. This is clear evidence that 
the calculated net charges on Cu and 0 are somewhat too small, and 
that the calculated electronic structure is more representative of 
the surface than of the bulk of the material . Higher binding 
energies for the 0 Is level are in fact observed in surface- 
sensitive experiments, which are in excellent agreement with the 
calculated value of 533 eV. For the different inequivalent Cu and 
0 atoms, no significant variations in the core-electron energies 
are found either in the experiment or in the calculation. 

The unit -cell superconductor calculations have been carried 
out on . the Cray supercomputers at the Numerical Aerodynamic 
Simulation (NAS) facility at Ames Research Center, whose support is 
gratefully acknowledged. 
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7 . COMPUTATIONAL METHODS 


Xa Methods . - The results for the low-work- function surfaces 
and high- temperature superconductors discussed in this report are 
based on the Xa scattered-wave approach. Fully relativistic (Case 
and Yang, 1980) and quasirelativistic (Cook and Case, 1989) 
computer programs have been employed for this research. In the 
quasirelativistic version, the relativistic mass-velocity and 
Darwin corrections to the energy are retained, but spin-orbit 
interaction is no longer taken into account. This reduces the 
computational effort considerably, and has enabled calculations for 
clusters with low symmetry and the study of adsorbate interactions 
as a function of coverage. 

The original Xa scattered-wave programs have been 
significantly improved. In the quasirelativistic program, the 
input /output has been changed to conform with the fully 
relativistic version. An interpolation scheme for structure- factor 
matrices (G) has been implemented which doubled the computational 
speed of both programs on the Cray supercomputers. Because the 
off-diagonal elements of the G matrices show generally a smooth 
variation with energy, they can be precalculated for a specified 
energy grid, stored on disk, and interpolated during the iterative 
process. In the final iterations calculated G matrices are used 
without interpolation. Another new useful feature allows one to 
automatically compute all unoccupied electronic states for a given 
system by employing an in-line energy search procedure. 

TABLE XXII, 

Net atomic charges q (in electrons) , dipole moments n, and dipole 
energies for BaO at R e = 3.666 a 0 = 1*940 A and Ba0/W 9 (100) with 
R(Ba-O) = 4.9 a 0 - 2.6 A and R (Ba-surface) = 6.0 a 0 -3.2 A . 


Property 

BaO 

quasi / fully rel. 

BaO/W 9 

quasi / fully rel. 

q (Ba) 

+ 0.82 

+0.85 

+1.50 

+ 1.51 

q (0) 

-0 . 82 

-0.85 

-0.15 

-0.24 

q (Wj) 



-0.08 

+0 . 02 

q(Wn) 



-0.00 

-0.33 

q(w XII ) 



-0.25 

-0.25 

fi (Debyes) 

9.5 

7.0 

38.4 

37.5 

E Ba (eV) 

4 . 0 

6.1 

6.1 

6 . 4 
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A comparison between quasi- and fully relativistic results for 
BaO and Ba0/W 9 (100) is made in Table XXII. The calculated charge 
distributions are in reasonable agreement, although there are some 
distinct differences in the net charges of 0, W x and W I: for BaO/W 9 . 
The dipole properties differ quite significantly for BaO, but they 
are similar for BaO/W 9 at the two levels of calculation. Because 
the dipole energy E Ba is defined as q 3 (Ba)//x, small variations in 
q (Ba) have a relatively large effect on the value of E^. 

Quasi- and fully relativistic results for BaO/M 4 W 1 W 4 clusters 
with M = W, Os, and Pt are compared in Table XXIII. It can be seen 
that the dipole moments at the quasirelativistic level (in 
parentheses) are slightly larger, by 2-5%, and that the dipole 
energies are 5-9% smaller than the corresponding fully relativistic 
results. The calculated data are otherwise very similar. 

The net charges for Ba and 0 on the different substrates vary 
only slightly (by 0.1 electrons or less), while the major charge 
redistributions occur inside the metal substrate and in the 
opposite direction of the Ba- substrate charge transfer. The 
consequences are a steady decrease in the dipole moments and an 
increase in the Ba binding energies in the series. 


TABLE XXIII. 

Comparison of fully and quasirelativistic (in parentheses) Xa 
results for BaO/M 4 W 1 W 4 clusters with M = W, Os, and Pt . R(Ba-O) = 
4 . 9 a 0 « 2 . 6 A. 


Property 

BaO/W 4 W x W 4 

Ba0/0s 4 W 1 W 4 

BaO/Pt^V^ 

q (Ba) 
q (0) 

q (Substrate) 

+ 1.51 
-0.24 
-1.27 

(+1.50) 

(-0.15) 

(-1.35) 

+ 1.58 
-0.17 
-1.41 

(+1.57) 

(-0.15) 

(-1.42) 

+ 1.58 
-0.14 
-1.44 

(+1.52) 

(-0.25) 

(-1.27) 

q (Mj) 

q(w„) 

q(W XII ) 

+ 0.01 
-0.33 
-0.25 

(-0.08) 

(-0.00) 

(-0.25) 

-0.21 

-0.09 

-0.12 

(-0.25) 

(+0.12) 

(-0.13) 

-0.29 
+ 0 . 07 
-0 . 08 

(-0.27) 

(+0.18) 

(-0.09) 

/ i (Debyes) 
Eb* (eV) 

37.5 

6.4 

(38.4) 
( 6.1) 

29.5 

9.2 

(31.1) 

J 8.6) 

26.6 

10,3. 

(25.8) 
( 9.4) 
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Benchmark calculations have been carried out on the different 
Cray supercomputers at NASA Lewis and Ames. In Table XXIV the 
timings are given for one iteration of the specified fully 
relativistic Xa calculations. The Cray-2 has only 60-70% of the 
speed of the Cray X-MP (in spite of a faster clock time, but due to 
a longer memory fetch time) , while the Cray Y-MP is faster than the 
X-MP by 30-40% for the present code. Without vectorization the CPU 
time increases on the Cray Y-MP by a factor of 2.8, and without 
optimization by a factor of 4.8. For the present applications, the 
Cray Y-MP is clearly the favored machine. 


TABLE XXIV. 


CPU timings on different Cray supercomputers for Xor calculations, 
in seconds, and relative speedup/down for the Cray- 2 and Y-MP 
relative to the Cray X-MP (in parentheses) . 


Calculation 

Cray Y-MP 

Cray X-MP 

Cray- 2 

BaO/W 9 

18.5" 

(1.0) 

23.6" (1.3) 

38.0" 

(2.1) 

BaO/W 9 + G-Mat. 

53.0" 

(1.0) 

73.5" (1.4) 

108.8" 

(2.1) 

Y 2 Ba 2 Cu 12 0 18 

190.0" 

(1.0) 

257.4" (1.4) 

390.1" 

(2.1) 


FLAPW Method . - The linearized augmented plane wave (LAPW) 
method is among the most accurate methods for performing band- 
structure calculations for solids or surfaces. Exchange and 
correlation are treated by using the local density approximation 
(LDA) to density-functional theory. Several forms of LDA 
potentials exist in the literature, the simplest being Slater's Xa 
potential. Relativistic effects can be included in a scalar- 
relativistic treatment for valence states and a fully relativistic 
one for core states. 


The LAPW method is a procedure for solving the Kohn-Sham LDA 
equation for the ground- state density, total energy, and 
eigenvalues (bands) of a periodic many-electron system by employing 
a basis set which is highly adapted to the problem. The adaption 
is achieved by dividing the unit -cell into atomic spheres, 
interstitial, and vacuum regions. Inside the atomic spheres linear 
combinations of radial functions times spherical harmonics Y lm (r) 
are used, and plane waves exp(ikr) in the other regions. The 
solutions to the Kohn-Sham equation are expanded in this basis and 
the coefficients are determined by the Rayleigh-Ritz variational 
method . 
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In its general form the LAPW method expands the potential as 
follows : 


E im v im( r ) Y lm (r) inside sphere, 

V (r) = 

E k V(k) exp(ikr) else. 

Thus no shape approximations are made, and the method is then 
called a full-potential (FLAPW) method. The muffin-tin 

approximation corresponds to retaining only the 1=0 and m=0 
components in the first expression for V(r) and only k=0 in the 
second. This procedure corresponds to taking the spherical average 
inside the spheres and the volume average in the other regions. 

A FLAPW program for solids (Blaha et al . , 1990) has been 

obtained from Prof. Trickey at the University of Florida in 
Gainesville. It has been implemented on the Cray Y-MP at Lewis and 
some tests have been run for body- centered- cubic tungsten. The 
program consists of several independent packages. After a starting 
potential has been generated, the self-consistency cycle involves 
the following steps and programs: 

1. Enter with density. 

2. LAPWO : Generate potential from density. 

3. LAPW1: Compute valence bands (eigenvalues and eigenvectors) . 

4 . LAPW2 : Compute valence density from eigenvectors . 

5. CORE : Compute core states and densities. 

6. Mixer: Mix input and output densities. 

7 . Loop to 2 . 


During his 1990 Summer Faculty Internship, Prof. K.S. Chu has 
developed a theoretical model for the treatment of surface 
relaxation effects within the FLAPW scheme by use of a deformation 
potential , which is based on elastic theory and lattice dynamics. 
This model and the implications for its implementation have been 
discussed with the Principal Investigator. 
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ft. CONCLUSIONS 


The quantum chemical cluster approach has been demonstrated to 
be a viable tool for the study of the electronic structure of low- 
work- function surfaces and high- temperature superconductors. 


Very reliable and consistent results have been obtained for 
the minimum work functions of different model cathode surfaces. 
The calculated Ba binding energies are in very good agreement with 
experimental energies of desorption. The optimum coverages depend 
strongly on the magnitude of the charge transfer at the surface, 
and there are indications that the calculated surface dipoles are 
too large and the resulting optimum coverages too small in some 
cases. The derived work functions, however, do not seem to be 
affected. 


For Ba and Cs on W(100) excellent results have been obtained 
for the minimum work functions and the adsorbate binding energies. 
In the case of Cs/W(100), the 0.25 monolayer coverage at the work 
function minimum is accurately reproduced. For Ba and BaO on 
W(100) , optimum coverages of -0.15 monolayers have been obtained, 
whereas higher coverages of up to 0.5 monolayers are expected. For 
BaO/W(110), on the other hand, the calculated optimum coverage of 
0.1 monolayer is consistent with experimental results. 


As a model for B-type cathode surfaces, it is found that the 
upriqht geometry for BaO on W(100) as suggested by Tuck and Norman 
et al. (1986 and 1987) provides a correct value for the B-type work 
function in relation to other low-work-function surfaces. The 
validity of this model has however been questioned by Shih et al. 
(1988) , who find a shorter Ba-0 distance for B-type surfaces and 
suqqest a bridge-bonded BaO structure. While a shorter Ba-0 
distance gives better results in comparison with photoemission 
data, the assumption of a bridge-bonded structure leads to higher 
work functions in the investigated cases. The calculated work 
function of the BaO/W(110) surface is also higher than as observed 
for B-type cathodes, so that this surface is not expected to 
contribute significantly to the emission. The correct BaO geometry 
and the specific surface that are responsible for the high emission 
from operating B-type cathodes remain, therefore, uncertain. 


The work function of BaO on hep Os (1010) is found to be lower 
bv 0.2 eV and the work function of BaO on fee Pt(lll) higher by 0.4 
eV than that for the model B-type Ba0/W(100) surface. These 
results are in excellent qualitative agreement with the observed 
emission enhancement for Os -coated cathode surfaces and t e 
degradation for Pt-coated surfaces. The variation m the work 
function is found to have its origin in the different substrate 
crystal structures and not the different electronegativities or 
work functions of the coating metals themselves. 
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•ff the calculated trends for the pure substrates can be 
verified also for M-type (alloy) surfaces, the credibility of the 
employed theoretical model would be significantly increased, and 
predictive capabilities could be achieved. Since the lowest work 
functions for dispenser cathodes are found for substrates with the 
hep crystal structure, the interplay between the magnitude of the 
surface dipole and the depolarization is expected to explain also 
the low work functions for other surfaces with the same crystal 
structure, like Os-W and Ir-W. 


The valence and core electronic structure of the 90K 
superconductor YBa 2 Cu 3 0 7 is found to be well described with atomic 
clusters of the size of the unit cell of the material . The 
experimentally observed intensity distribution and width of the Cu- 
0 valence band are correctly described. The valence band observed 
in photoemission spectra is shifted, however, by about 1.5 eV to 
higher binding energies. This effect appears to be present in all 
density- functional calculations. Final-state energy corrections in 
the localized Cu 3d states are found to contribute 10-20% to the 
observed shift, while true many-electron effects have to made 
responsible for the majority of the shift. The calculated core- 
electron binding energies for Cu, O, and the 3d states of Y and Ba 
are, after inclusion of the final-state corrections, in excellent 
agreement with experiment. The charge distribution, as obtained 
with the cluster approach, is believed to be representative of the 
surface of the superconducting material. 
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NAS Technical Summary, March 1988 - February 1989 


ELECTRONIC STRUCTURE OF SUPERCONDUCTORS 

Wolfgang Mueller, Principal Investigator 
Analaiom, IncVNASA Lewis Research Center 

RESEARCH OBJECTIVE: 

*n» primary objective of this work is to establish the validity of the quantum chemical cluster approach for the description of the electronic structure of the new 
high 'temperature superconductors. This will justify theoretical analyses* enable predictions for materials properties* and contribute to a better undemanding of 
the mechanism responsible for high -tempera turn superconductivity. 

APPROACH: 

Large cluster calculations are performed using the self -consistent relativistic scattered-wave method* which is based on the Dirac wave equation. The method 
employs the multiple scattering technique and Slater's statistical approximation for electron exchange and correlation. 

ACCOMPLISHMENT DESCRIPTION: 

The fust fully relativistic unit-cell cluster results have been obtained for YBa 2 Cu 3 O 7 . The accompanying figure shows the density of states (DOS) for a 
XfcB^Cu jjOi g duster. Tbe DOS is generated from the calculated energy levels which are indicated as small vertical lines on the horizontal axis. The peaked 
distributions were obtained by broadening the cluster energy levels with Gaussians using width parameters of 0.05 and 0.50 eV, respec- 
tively. The calculated DOS is in good agreement with solid-state results and experimental photoemission spectra. The calculated net atomic charges in 
YBa 2 Cu 3 07 are, approximately* for Y: +U; Ba: +0.6; Cul: +0.6; Cu2: +0.8; 01, 02, 03: -0.7; and 04: -O. 6 . A typical calculation requires about 5 Cray-2 
hours and 2 megawords of memory. 

SIGNIFICANCE: 

A thorough understanding of the electronic structure of the existing high -temperature superconductors will be crucial for the development and discovery of new 
materials with higher critical temperatures and current densities. The cluster approach is particularly well suited to provide direct insight into important molecu- 
lar parameters through the analysis of the orbital structure at the Fermi energy. 

FUTURE PLANS: 

The code will be further used to investigate valence and core-level properties, relativistic electronic effects* and the coupling between Cu-O planes in other 
high-temperature superconductors. Based on the results of the scattered-wave study, calculations will be designed which employ more elaborate ab initio 
techniques. 

PUBLICATIONS: 

i Relativistic Ouster Results for YBa 2 Cu 3 <> 7 , Wolfgang Mueller. To be presented at the International Conference on Materials and Mechanisms of Super- 
conductivity: High-Temperature Superconductors, July 23-28* 1989* Stanford, California* and to be published. 
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NAS Technical Summary, March 1989 - February 1990 


Electronic Structure of Superconductors 


Wolfgang Mueller, Principal Investigator 
Anaiatom, Inc ^ ASA Lewis Research Center 


Vteeearch Objective 

To investigate the electronic structure of high-temperature 
superconductors with the quantum chemical duster approach, 
perform theoretical analyses for critical materials properties, 
and contribute to a better understanding of the mechanism 
responsible for high-temperature superconductivity. 
Approach 

Large-duster calculations are carried out using the self- 
consistent, relativistic scattered-wave method, which is based 
on the Dirac wave equation. The method uses the multiple 
scattering technique and Slater's statistical approximation for 
electron exchange and correlation. 

Accomplishment Description 
Putty relativistic calculations have been performed for 
Y-8a-Cu-0 using two different clusters: Yj8a2Cui20i8 with 
three Cu-0 layers, and YBajCuneOa with four Co-0 layers. 
The accompanying figure shows the valence and conduction 
electronic structure per cell as derived from the 

two dusters, and a decomposition Into atomic contributions. 
The prominent double peak just below the Fermi level, and 
other features that are observed in photoemission experi- 
ments, are clearly resolved. For YjBa^uijOie, final-state 
photoemission relaxation contributions have been calculated 


for about 25 valence and core states. Final-state SffecHt * 
(relative to the duster Fermi level) are found to be smalt within 
the valence band, but significantly affect the core levels, 
where the remaining errors are 1 to 5 eV. A typical calculation 
requires about 5 Cray-2 hours and 2 megawords of memory. 
Significance 

The investigation of spin-orbit, relaxation, and correlation 
effects in duster calculations Is important for understanding 
the discrepancies between band structure results and 
photoemission data. Large clusters are shown to provide^ 
good description of the electronic structure of these materials. 


Future Flans 

Extended duster and plane-wave calculations wlflbapar- 
(onwd for tNn film# of Y-B*-Co-0 on various substrates. 


puMieaimw — — > — 

1 . Mueller. W. "Relativistic Cluster Results ter 
Physic* C 162-164 (1660): 1357. 

2. Mueller, W. *Valenc# and Cora Electronic Structure of the 
90K Superconductor Y-Ba-Cu-O.* Presented at the 10th 
Canadian Symposium on Theoretical Chemistry, Banff, 
Alberta. Aug. 1969. 


DENSITY OF STATES FOR Y 2 Ba 2 Cu B 0« 



DENSITY OF STATES FOR YBa 2 Cu«0» 



Density of states for (top) YjBazCu-ijO^ and (bottom) 
YBa2Cu 16 02o 



Atomic contributions to the density of states. 
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NAS Technical Summary- 


March 1990 - February 1991 


Electronic Structure of Superconductors 

Wolfgang Mueller, Principal Investigator 
Analatom, IncJNASA Lewis Research Center 


Rtsurch Objective 

To Investigate the electronic structure of high-temperature 
superconductors with the quantum chemical-cluster approach, 
to perform theoretical analyses for critical materials properties, 
and to contribute to a better understanding of the mechanism 
responsible for high-temperature superconductivity. 


Approach 

Large-cluster calculations are carried out using self-consistent 
ab initio and scattered-wave methods, which are based on the 
Schroedinger and Dirac wave equations. 

Accomplishment Description 
A comparison has been made among the results of fully, 
quasi-, and nonrelativistic scattered-wave calculations for 
Y-Ba-Cu-O; a representation of a Y^Ba^Cu^O^ cluster with 
three Cu-0 layers was used. The cluster is shown schemati- 
cally in the figures. The two figures depict the valence and 
conduction electronic structure per YBaXugOj cell as derived 
from the fully and quasl-relativistic calculations, respectively. 
The densities of states have been generated from the calcu- 
lated energy levels (which are indicated as small vertical lines 
on the horizontal axis) by Gaussian broadening with width 
parameters of 0.1 and 0.5 eV. Although the general electronic 
structure is similar in both calculations, especially with respect 
to the prominent double peak just below the Fermi level, there 
are significant differences in the shapes of the Individual 


features and their positions relative to the cluster Fermi level. 
The density of states from the quasl-relativistic calculation is in 
good agreement with quaskelativistic band structure results, 
which demonstrates the validity of the cluster approach. A 
typical calculation requires 5 to 10 Cray hours and 4 mega- 
words of central memory. 

Significance 

The investigation of relativistic, relaxation, and correlation 
effects in cluster calculations is important for a detailed 
understanding of the electronic structure and the coupling 
mechanism in high-temperature superconductors. Fully 
relativistic and correlated wave functions are required for high 
accuracy. 

Future Plana 

Embedded cluster and plane-wave calculations will be carried 
out to investigate differences in the surface, bulk, and inter- 
face electronic structure of Cu-0 superconductors and films 
on various substrates. Electronic coupling parameters will be 
determined using the ab initio approach. 

Publications 

Mueller, W. *A Comparison of Theoretical and Experimental 
Photoemission Data for YBajCuoO^’ Presented at the 7th 
North Coast Symposium of the American Vacuum Society, 
Cleveland, OH, May 1990. 




(a) Fully relativistic. 


(b) Quasl-relativistic. 


Electronic structure of Y-Ba-Cu-O. 
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RELATIVISTIC CLUSTER RESULTS FOR YBa 2 Cu 3 0 7 * 


Wolfgang MUELLER 

Analatora, Inc., 29904 Sycamore Oval, Westlake, Ohio 44145, USA 


The valence density of states, charge distribution, and core-level binding energies for a unit-cell 
cluster of Y-Ba-Cu-0 are compared with solid state and experimental results. The cluster is shown to 
be a reasonable model, in particular for surface properties. Cu 3d spin-orbit effects and final state 
contributions in the photoeroission process are Investigated for the first time. 


1. INTRODUCTION 

It is known 1 ' 2 that density-functional solid 
state calculations do not provide an accurate 
description of the Cu 3d - 0 2p derived density 
of states (DOS) in the copper oxide super- 
conductors. Although the general structure of the 
photoemission spectrum Is well reproduced, the 
experimental spectra are shifted by 1-2 eV to 
higher binding energy, and the DOS at the Fermi 
level is much smaller than calculated. This dis- 
crepancy is attributed to many-electron effects 
in the highly localized Cu 3d states, which can 
be treated, in principle, with high accuracy in 
cluster calculations. As a test for the validity 
of the cluster approach, the electronic structure 
of large clusters is investigated, still within 
the density functional approach, but with 
consideration of Cu 3d spin-orbit effects and 
photoemission final state contributions. In these 
respects, the present calculations go beyond 
standard solid state techniques. 

2. METHOD 

Large cluster calculations are performed for 
YBa2Cu 3 0 7 using the self-consistent, fully 
relativistic Xa scattered-wave approach. 3 Final 
state contributions are calculated by employing 
Slater's transition state method. The results 
presented here are based on a Y2Ba2Cui 2 0i8 

* Work supported by NASA Lewis Research Center, 
Cleveland, and the NAS program at NASA Ames. 


unit-cell cluster, which consists of a central 
Cu 4 (l)0 2 (l) plane, with Ba0 4 (4), Cu 4 (2)0 4 (2,3) , 
and Y planes on both sides, t r 


3. VALENCE DENSITY OF STATES 

The total and partial densities of states for 
YBa2Cu30 7 as derived from the 34-atom cluster are 
shown in Fig. 1. The peaked and smooth 
distributions for the total DOS per cell were 
obtained by broadening the Indicated cluster 
energy levels with Caussians using width 



el CURE 1. Total (per YBa 2 Cu 3 0 7 cell) and partial 
(atomic) densities of states as derived from a 
Y2Ba 2 Cui20i8 cluster model for Y-Ba-Cu-0. 


0921-4534/89/S03.50 © Elsevier Science Publishers B.V. 
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parameters of 0.05 and 0.5 eV, respectively. The 
calculated DOS is in reasonable agreement with 
solid state results.* It shows the two major 
maxima observed in photoemission spectra, 
separated by about 2 eV, and additional features 
at higher and lower binding energies. The cluster 
DOS is shifted, however, by about 1.5 eV to lower 
binding energy, very similar to the DOS from 
solid state calculations. Apart from this shift, 
the partial DOS from the cluster show major 
copper contributions near the center of the 
double maximum and oxygen contributions more to 
Its sides, in good agreement with experiment. 1 The 
Cu(l,2) 3d states near -2 eV are strongly 
localized and clearly show the spin-orbit 
splitting of about 0.3 eV. The Cu(l)-0(1 ,4) and 
Cu(2)-0(2,3) states are otherwise significantly 
hybridized. Without spin-orbit interaction, the 
Cu-derived DOS from solid state calculations 
exhibits much sharper peaks. 

The calculated direct energy shift due to the 
relaxation of the charge distribution upon 
ionization, which is the final state contribution 
in the photoemission process, is very similar 
Jl.4 eV) at the lower and upper ends of the 

valence band of the cluster. What counts in 
comparison with photoemission spectra is the 

shift of different states relative to the 
calculated shift at the Fermi level. The presence 
of an almost Identical shift at both ends of the 
valence band therefore indicates that its width 
is not affected by final state contributions. The 
largest relative shifts of up to 0.2 eV are 

obtained for the strongly localized Cu 3d states 
near -2 eV. Although quite small, these final 
state effects account for about 10% of the 

observed shift of the experimental double peak 
relative to the calculation. 

4. CHARGE DISTRIBUTION AND CORE-LEVEL ENERGIES 

The analysis of the charge distribution shows 

a Cu 3d population of 9.3 for Cu(l) and 9.6 for 

Cu(2). These values are in very good agreement 


with the estimate of 9.45 from the Larsson- 
Sawatzky model. 4 5 Because of the presence of a 
strongly covalent electronic structure and the 
division of charge among overlapping spheres and 
from the interstitial space, the Cu atoms 
formally acquire also substantial amounts of 
valence s and p charge, so that their net charges 
become +0.6 and +0.8 electrons, respectively. 
With this subdivision of charge, the excess 
charge on oxygen, -0.7 for 0(1, 2, 3) and -0.6 for 
0(4), as well as the ionicity of Y (+1.1) and 
Ba(+0.6) are small. 

A comparison of calculated core-level binding 
energies with experiment is therefore of 
particular importance. The calculated Cu 2p 1/z 
energies, for both sites and Including the final 
state contribution, are within 952*1 eV (expt. 
953 eV), 6 and Cu 2p 3 / 2 within 932*1 (933 eV). 6 - 7 
The Y and Ba 3d energies are also within 1-2 eV 
of experimental values. The calculated 0 Is 
energies for all four sites are within 533*1 eV. 
Peaks are observed at 533, 531, 6 and 529 eV; 7 
however, only the latter peak is attributed to 
the bulk material. The fact that the cluster 0 Is 
energies are all In the vicinity of the higher, 
surface-sensitive values indicates that the net 
charge on 0 in the isolated cluster is, due to 
the missing Madelung field, too small. 

In conclusion, the electronic structure in Y- 
Ba-Cu-0 is both ionic and covalent in character. 
The cluster valence DOS is in good overall 
agreement with solid state and experimental 
results. The charge distribution and core-level 
energies Indicate, however, that an Isolated 
cluster should be considered more as a model for 
the surface of the material than for the bulk. 


REFERENCES 

1. J.-M.Mariot et al., Z.Phys.B 75,1(1989). 
J.A.Yarmoff et al., Phys. Rev. ff~36, 3967(1987) . 

?,*5o£? Se and c - Y * Yan 9. J-Chem.Pliys. 72,3443 
(1980) • 

S.MassIdda et al., Phys.Lett.A 122,198(1987). 
A. Fujimori et al., Phys.Rev.B 3575814(1987). 
D.D.Sarma et al., Phys.Rev.B 3572371(1987). 

P. Steiner et al., Z.Phys.B 67719(1987). 


2 . 

3. 

4. 

5. 

6 . 
7. 


77 


..Jill 





pmmum pagi blank not filmed 

180 IEEE TRANSACTIONS ON ELECTRON DEVICES, VOL. 36, NO. I, JANUARY 1989 

Electronic Structure of BaO/W Cathode Surfaces 

WOLFGANG MULLER 


Abstract— The local electronic structure of the emissive layer of bar- 
ium dispenser thermionic cathodes is investigated theoretically using 

— the relativistic scattered-wave approach. The interaction of Ba and O 
with W, Os, and W-Os alloy surfaces is studied with atomic clusters 
modeling different adsorption environments representative of B- and 

* s M-type cathodes. Barium is found to be strongly oxidized, while oxy- 

— gen and the metal substrate are in a reduced chemical state. The pres- 
ence of oxygen enhances the surface dipole and Ba binding energy rel- 

^ : ative to Ba on W. Model results for W-Os alloy substrates show only 

J j relatively small changes in Ba and O for identical geometries, but very 
large charge redistributions inside the substrate, which are attributed 
to the electronegativity difference between Os and W. If Os is present 
t ; in the surface layer, the charge transfer from Ba to the substrate and 
l ^ the Ba binding energy increase relative to W. The improved electron 

— emission from modern alloy surfaces could therefore possibly be ex- 
plained by higher coverages due to increased Ba binding energies. The 
different emission enhancement for different alloy substrates is likely 

■ to be related to a competition between electronegativity effects and the 
availability of unoccupied rf-type orbitals in the substrate. 

I. Introduction 

m T^HERMIONIC cathodes have many important appli- 
X cations in the vacuum tube industry. In the field of 
1 4 power microwave devices [1], cathodes are major com- 
ponents in klystrons and traveling-wave tubes used as ra- 
dio frequency signal amplifiers for communications. Other 
applications include high-resolution radar systems, plasma 
heating for fusion research, and free-electron lasers. 
z — The basic and still most widely used type is the barium 

dispenser cathode, originally designed in 1955 by Levi 
[2] in which barium calcium aluminate is impregnated into 
= a porous tungsten matrix. These so-called B-type cath- 

odes have a work fifnction, of about 2 eV and operate at 
around 1400 K. The search for a cathode with a lower 
operating temperature and'resulting longer life led in the 
1960’s to the development of M-type cathodes [3] in 
which the tungsten surface is coated with a film of rhen- 
ium, osmium, osmium-ruthenium, or iridium. A recently 

— improved reservoir cathode developed by NASA and Var- 
ian [4], with a W-Os alloy emitting surface, has the po- 
tential to operate at current densities of 1-100 A/cm 2 and 

_ is expected to have a very long life. Contemplated appli- 
cations are high-frequency/high-power tubes for space 
communications and microwave power transmission from 
satellites. 

Extensive experimental efforts using modem surface 
analytical techniques like X-ray photoelectron (XPS) and 

Manuscript received July 15, 1988; revised September 29, 1988. This 
1 work was supported by the NASA Lewis Research Center under Contracts 
NAS3-24744 and NAS3-25085. 

The author is with Analatom, Incorporated, Westlake, OH 44145. 

IEEE Log Number 8825212. 


Auger electron spectroscopies (AES) have been under- 
taken to characterize the electronic structure of the cath- 
ode emissive layer. The adsorbate structure is roughly de- 
scribed as a monolayer of barium above oxygen on the 
tungsten(-alloy) substrate (see, e.g., [5]). The detailed 
chemical and electronic structure of the cathode surface, 
however, has been found to be complex and is experi- 
mentally not well understood. Recent surface-extended X- 
ray absorption fine structure (SEXAFS) measurements [6], 
[7j of real cathode surfaces have resulted in conflicting 
assignments of atomic distances and even the geometrical 
arrangement of the Ba and O atoms in the emissive layer 
of impregnated tungsten cathodes. 

Although the improved electron emission from modem 
low-work-function, cathode surfaces is experimentally 
well documented, basic structural and electronic surface 
properties and, in particular, and mechanism for the emis- 
sion enhancement in alloy cathodes are not understood 
[8]. Because of the ever increasing demand for cathodes 
with higher current densities and longer life, a thorough 
theoretical understanding of the electronic properties of 
the surface at the atomic level is very much needed. 

The work reported here represents the first quantum 
chemical investigation, based on the cluster model ap- 
proach, of the emissive layer of barium dispenser 
thermionic cathodes. A preliminary account of the major 
results obtained from a BaO/W 9 cluster has been given 
elsewhere [9]. 

The paper is organized as follows: The relativistic clus- 
ter model approach and computational details aje briefly 
discussed 'in Section II. The calculated electronic orbital 
structure of BaO/W ( 100) is compared with experimental 
photoemission data in* Section III. In Section IV, the 
charge distribution at the surface is analyzed and the re- 
sulting dipole properties are investigated. In Section V the 
first computational results for models of W-Os alloy sur- 
faces are presented. The conclusions from the present the- 
oretical study are summarized in Section VI. 

II. Relativistic Cluster Model Approach 

Electronic structure calculations are performed for 
atomic clusters representative of the BaO/W(-Os) sur- 
face complex by employing the relativistic Xa scattered- 
wave method. 

The fully relativistic self-consistent field, Xa scattered- 
wave method is based on the one-electron Dirac equation 
and has been developed by Yang and Case [ 10]— [ 12] . The 
general formulation follows the nonrelativistic approach 
of Slater and Johnson [13], [14] but includes all relativ- 


00 1 8-9383/89/ 0100-01 80$0 1 . 00 © 1989 IEEE 



MULLER: ELECTRONIC STRUCTURE OF BaO/W CATHODE SURFACES 

istic electronic effects: the mass velocity and Darwin cor- 
rections to the energy, and spin-orbit interaction through 
the use of four-component electronic wave functions in 
the framework of double group theory. Electron exchange 
and correlation effects are taken into account by employ- 
ing Slater’s Xa potential [15]. The program uses overlap- 
ping atomic spheres and thus goes beyond the basic muf- 
fin-tin approximation for the potential. The method 
contains no adjustable molecular parameters. 

The substrate is modeled with atomic clusters repre- 
senting the (100) face of the body-centered-cubic (bcc) 
crystal lattice of tungsten. Although real cathode surfaces 
are polycrystalline and ( 100) may not give the dominant 
contribution to the overall cathode emission, this surface 
has been chosen because of its high symmetry, which 
considerably facilitates the calculations. In a model for 
BaO adsorption on osmium, the (1010) face of the hex- 
agonal-close-packed (hep) lattice has been chosen. W-W 
and Os-Os distances as present in the bulk are used, with 
lattice parameters of 3.165 A for W and 2.74 and 4.32 
A for Os [16]. 

Atomic clusters with between 5 and 14 substrate atoms 
are investigated. The larger clusters with 9 and 14 sub- 
strate atoms are shown in Fig. 1. The smallest cluster, 
W 5 , is derived from the body-centered-cubic W 9 cluster 
by considering only the 4 symmetry-equivalent atoms in 
the first layer and their nearest neighbor in the second 
layer. W 10 contains in addition to W 9 the nearest neighbor 
in the fourth layer. In W 14 two W 9 units are combined. 
Mixed W-Os substrate clusters are derived from W 9 by 
selectively replacing W atoms by Os. 

For the bcc (100) clusters, oxygen adsorption is as- 
sumed in the fourfold hollow site with barium on top of 
oxygen. In the hep Ba0 2 /Os 9 cluster, barium occupies 
the hollow site and is bridge-bonded to two oxygen at- 
oms. A Ba-O bond length of 4.9 Bohrs ( * 2.6 A) is used 
in all clusters, as determined by Norman et al. [6] from 
SEXAFS measurements for BaO on W and W-Os cathode 
surfaces. (Note, however, that Shih et al. [7] in a recent 
paper find 2.3 A for the W substrate). The Ba-W distance 
used is 3.9 A. This value is based on a fit of surface 
models to the Norman data and has been reported by Tuck 
et al. [17], The Ba-Os distance used in Ba0 2 /0s 9 is 3.75 
A and taken the same as determined for BaO/ W-Os [6]. 
This leads to the following heights above the surface (in 
Bohrs = 0.5292 A): and 1.1 for O 6.0 for Ba in all bcc 
clusters, and 5.2 for Ba and 1.5 for O in Ba0 2 /0s 9 . 

Generally, the Ba electrons up to 4 d, the O Is electrons, 
and the W and Os electrons up to 5 p and including 4/are 
considered core electrons, for which spherical atomic po- 
tentials are assumed. The Ba 5s and 5 p semi-core elec- 
trons and the remaining valence electrons are treated in 
the full cluster potential. The Ba 2 0 2 /W !4 , Ba0 2 /Os 9 , and 
BaO / W 2 Os 2 W | Os 2 W 2 clusters have C 2v point group sym- 
metry, all other clusters C 4v symmetry. The radii of the 
atomic spheres are chosen as 88 percent of the sphere radii 
determined from the superposition of nonrelativistic 
atomic charge densities [18]. Xa values of 0.74 for O and 
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Fig. 1. Cluster models for BaO/W(Os). Mixed W-Os substrate clusters 
are derived from BaO/W y by selectively replacing W by Os. 

0.70 for Ba, W, Os, and the outer sphere are used. For 
the largest cluster, Ba 2 0 2 /W )4 , a total of 1164 electrons 
are treated self-consistently, which leads to 110 different 
core states and 58 valence states. With partial waves in- 
cluded up to l = 1 for O and / = 2 otherwise, the complex 
secular matrix for the semi-core and valence states of 
Ba 2 0 2 /W| 4 has a dimension of 161. 

III. BaO/W Electronic Structure 

Fig. 2 displays the calculated outer electron density of 
BaO/W 9 in a plane perpendicular to the surface along the 
[110] direction, as indicated in the atomic cluster. Con- 
tributions from the W 5 d, 6s, 6 p, O 2s, 2 p, and Ba 5s, 
5p, 5 d, 6s, 6 p electronic levels are included. The elec- 
tronic configuration of tungsten in the clusters is approx- 
imately 5d*6s l , including some 6 p character. The d or- 
bital structure on W with a characteristic density of four 
lobes is clearly visible in the figure. 

The electron density around the barium atom is largely 
spherically symmetric due to a high level of oxidation; 
what is seen is essentially the completely filled 5s and 5 p 
semi-core shell. The oxygen atom is embedded in the first 
tungsten layer and exhibits a large interaction in particular 
with the second-layer tungsten atom. The population 
analysis shows that the oxygen orbitals are strongly mixed 
with tungsten orbitals: The O 2s orbital has 10 percent 
and the O 2 p orbitals have each about 30 percent of W 
character, mainly in W 5 d orbitals. 

The orbital energy structure of BaO/W 9 gives rise to 
the density of states (DOS) shown in Fig. 3, in which the 
discrete cluster levels have been broadened by a Gaussian 
distribution with a half width of 0.3 eV . The arrows in- 
dicate observed energies from the photoemission study of 
Mueller et al. [19] for a BaO overlayer on W( 100). The 
calculated density of states below the Fermi level down 
to -6 eV originates from electronic states related to oc- 
cupied W 5 d and 6s orbitals. The BaO/W 9 spectrum 
shows the well-known W surface states at the Fermi level, 
the prominent surface resonance state at -1.3 eV (ex- 
perimentally at — 1.9 eV), and a three-peak structure be- 
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Fig. 2. Cluster model and calculated outer electron density of BaO/W 9 in 
a plane perpendicular to the surface. 



Fig. 3. Calculated density of states for BaO/W 9 . The vertical lines on the 
horizontal axis indicate the discrete electronic levels of the cluster and 
the arrows experimental photoemission energies for BaO/W(100) [19]. 


tween -3 and -6 eV. The surface resonance state and 
the two peaks between -3 and -4 eV are shifted relative 
to clean W 9 by 0. 2-0.3 eV to higher binding energies, 
while the peak below -5 eV is shifted by 0.4 eV to lower 
binding energy. The intensity of the peak near -4 eV is 
increased relative to clean W 9 . The three-peak structure 
is weakly seen in the experimental spectrum for clean 
W ( 100) but is not well resolved for BaO/W ( 100) f 19]. 
The intensity increase around —4 eV, however, is clearly 
observed. In the experimental spectrum, the structure at 
higher energies overlaps significantly with a broad O 2 p 
band with a maximum at -5.9 eV, which is shifted to 
— 8.5 eV in the theoretical spectrum. The calculated peak 


positions are otherwise in good agreement with experi- 
ment (within 0. 5-1.0 eV). 

The calculated electronic energies for three different 
clusters are compared in Table I with the experimental 
photoemission energies from Mueller et al. [19] and cal- 
culated energies based on the full-potential linearized 
augmented plane wave (FLAPW) thin film method from 
Hemstreet et al. [20]. The model used in the FLAPW cal- 
culation was a c(2 X 2) BaO overlayer on W( 100) with 
fl(Ba-O) = 2.375 A and R(Ba-W) = 3.617 A. 

A close similarity is observed for the calculated elec- 
tronic energies, which differ for the various orbitals in the 
quite different clusters by no more than 0.5 eV. This re- 
sult is particularly significant, since the energies in Table 
I are binding energies relative to the highest occupied 
(cluster Fermi) level, which changes from 4.8 eV for 
BaO/W 5 to 5.4 eV for BaO/W 9 and to 7.3 eV for 
Ba 2 0 2 /W 14 , i.e., by as much as 2.5 eV. Cluster size ef- 
fects thus appear to be not very crucial for electronic bind- 
ing energies. 

The comparison of the cluster results with experiment 
for binding energies less than 50 eV shows good agree- 
ment (within 1 eV) for W 5p 3 / 2 , Ba 5 p\/ 2 and 5p 3 / 2 , and 
O 2s . The Ba 5s level overlaps with W 4 / and no exper- 
imental value was given. The W 4 /electrons are more 
weakly bound in the calculation by about 3.5 eV, which 
could have resulted from the constraints imposed by the 
treatment as core orbitals. The calculated spin-orbit split- 
tings, however, of the W 4/ as well as the Ba 5 p levels 
are in excellent agreement with experiment, within 0. 1 eV 
of the total spin-orbit separations of 2 eV. 

The most significant discrepancy is in the O 2p binding 
energy, which is calculated about 2.5 eV higher than ob- 
served. This suggests that the height of oxygen above the 
W( 100) surface of 0.6 A as derived from the geometry 
proposed by Tuck et al. [17] may not be correct. For 
larger oxygen-surface separations a larger charge transfer 
to O and up to 2-3 eV smaller O 2 p binding energies are 
obtained [21]. The FLAPW calculation for a 0.1 A 
smaller oxygen-surface separation gives a split O 2 p peak 
with maxima at 7.8 and 8.8 eV, which is otherwise in the 
same energy range as the present results. The deviations 
of the FLAPW binding energies from experiment for the 
non spin-orbit split Ba 5 p levels and the O 2s level are 
3-4 times larger than for the cluster results. 

It should be noted at this point that the quoted orbital 
binding energies are, per definition, not equivalent to ob- 
served photoemission energies, but have to be corrected 
for self-energies and final state effects. These corrections 
are generally small for energies close to the Fermi level 
and are known to increase the calculated values for higher 
binding energies. The fact that the calculated Ba 5 p and 
O 2s, 2 p energies are higher than the experimental values 
for all clusters therefore indicates that the Ba + character 
is slightly too. high and the 0“ character too low. This is 
to be expected because the surface coverage in the largest 
cluster (Ba 2 0 2 /W 14 , which gives the best agreement with 
experiment) is only 1 /3 and in the photoemission exper- 
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TABLE I 

Electronic Enf.rgif.s Relative to the Highest Occupied (Fermi) Level 
for Selected Orbitals of BaO/W ( 100) in Electronvolts 
(The W energies given are those for the surface layer.) 


Orbital 

BaO/Vs 

BaO/U, 

Ba^O^/U x A 

Expt.[I9] 

FLAPU [201 

U 5p,, 3 

45.5 

45.8 

45.7 




35.0 

36.0 

36.0 

36.0 


B/7 

29. B 

30.0 

29.0 

33. 5 


4f 

27.5 

27.7 

27.6 

31.3 


Ba 5s 

31 . 0 

30.9 

30.5 



5Puz 

17. 5 

17. 3 

17. 1 

16.6 

13.0 


15.4 

15.2 

14.9 

14.6 

13.0 

0 2s 

22.9 

22.7 

22.5 

21 . 9 

25.0 

2p 

B. 7 

8.5 

0.3 


e.e p 7.8 


iment presumably 1 /2. At higher coverages, the net pos- 
itive charge on Ba will decrease because of dipole-dipole 
depolarization effects, and the net negative charge on O 
may increase because the oxygen-surface distance could 
be larger than assumed. There is also the possibility that 
the vertical arrangement for BaO on W( 100) is not cor- 
rect. 

IV. Charge Distribution and Dipole Properties 

In this section the charge distribution in BaO/W clus- 
ters will be analyzed. The results will be studied as a 
function of cluster size and compared with Ba/W in order 
to analyze the roles of the metal adsorbate and oxygen. 

The first basic quantity that can be calculated from the 
charge distribution is the surface dipole moment. Al- 
though the surface coverage in the clusters is in most cases 
equal to 1/4, the calculated dipole moments are not to be 
understood as the effective dipoles that determine the 
value of the work function for this coverage. Rather, the 
calculated dipole moments reflect initial dipoles, govern- 
ing the initial change of the work function with coverage 
in the low coverage regime. This is so because the inter- 
action between the dipoles from different adsorbates is not 
present in most of the clusters considered. The depolari- 
zation resulting from the repulsive interaction of two BaO 
adsorbate dipoles in Ba 2 0 2 /W| 4 is found to be very 
strong. The present results, therefore, emphasize funda- 
mental properties of an individual adsorbed molecule 
rather than properties of a more complete overlayer. The 
cluster results thus complement thin film calculations in 
which dipole properties and data for low coverages are 
not easily accessible. Symmetry considerations in thin film 
(slab geometry) calculations require the adsorbate to be 
placed on both sides of the film. The dipole moment for 
the system is, therefore, by definition equal to zero. 

A second basic quantity that can be derived from the 
charge distribution and calculated dipole moment is the 
dipole energy that is an estimate for the binding energy 
of an ionic adsorbate. In the classical picture, the net 
charge on the adsorbate interacts with its image inside the 


substrate to form an ionic adsorbate-substrate bond. The 
dipole binding energies obtained from the calculated ac- 
tual charge distributions will be compared with those ob- 
tained from the classical image charge concept and with 
experimental energies for desorption. 

The calculated properties for the different clusters are 
presented in Table II. These include the charge distribu- 
tion, the Ba 5 d orbital occupation, the dipole moment, 
and estimates for the Ba binding energy based on the elec- 
trostatic dipole and image charge concepts. Results for a 
stretched BaO molecule and for Ba on W are given for 
comparison, with the same distances used as in the 
BaO/W surface clusters. The roman subscripts on W de- 
note the different substrate layers. 

The charge distribution in the BaO/W clusters shows 
that electrons have been transferred from Ba to O and W . 
The charge transfer increases significantly from the two- 
layer W 5 to the three-layer W 9 cluster, but stabilizes very 
well with the four-layer W 10 cluster for the single-BaO- 
adsorbate systems. The differences in the individual net 
atomic charges are around 0. 1 electrons, with the largest 
variation in the substrate atoms. In the larger clusters, the 
first W layer is almost neutral, while the electronic charge 
from Ba has been transferred to the second substrate layer 
and to a large extent also to the third layer. Because the 
W 5 cluster does not include third-layer atoms, the charge 
distribution in the first two layers is completely different. 
The net charge on Ba in stretched BaO, Ba/W 9 , and 
BaO/W 5 is +1 electron. 

For the larger single-BaO-adsorbate clusters, the oxi- 
dation level of Ba is about +1.5. The Ba 5 d occupation 
is of the order of 0.3 electrons and the Ba 6 sp occupation 
about one half of the d occupation. This indicates a strong 
Ba 5 d contribution to the bonding. The net charge on ox- 
ygen is -0.2 electrons and surprisingly small. In the 
Ba 2 0 2 /W| 4 cluster with two BaO adsorbates in neighbor- 
ing sites, the net positive charge on Ba decreases to +1.2 
electrons and the net negative charge on O increases to 
-0.35 electrons. For larger oxygen-surface distances, for 
which the O 2 p binding energies would be in better agree- 
ment with experiment, the net charge on O may increase 
up to -0.5 electrons [211. 

The important general results from the analyses of the 
cluster charge distributions are summarized as follows: 
First, electronic charge is removed from Ba and found in 
form of a screened image charge inside the substrate clus- 
ter. Second, the net negative charge on O is small. And 
third, the tungsten substrate is not oxidized, it is in a re- 
duced metallic state with the surface layer of W atoms 
almost neutral and metal-like. The electronic structure of 
the emissive layer, therefore, is described as largely ionic 
between Ba and the O/W substrate, while oxygen forms 
covalent donor bonds to tungsten. The Ba-0 overlayer is 
completely different from either molecular or bulk barium 
oxide. According to the present results, the chemical 
structure is best characterized as a Ba salt of the O/W 
substrate. 

The large charge transfer from Ba into the substrate re- 
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TABLE II 

Net Charges <?( Electrons), Dipole Moments and Dipole and 
Image Charge Binding Energies E for Barium in Ba, BaO/W 
Clusters (R(Ba-O) = 4.9 Bohrs ^ 2.6 A, 
R(Ba-surface) = 6.0 Bohrs ~ 3.2 A.) 


Proper ty 

BaO 

Ba/V„ 

Ba0/U o 

BaO/U* 

BaO/W, o 

Ba 3 0 3 /U,« 

q ( Ba , tot ) 

♦ 1.03 

♦ 1.04 

♦ 0.90 

♦ 1.51 

♦1.47 

♦1.22*2 

q < Ba , 5d ) 

0. 76 

0. 53 

0. 76 

0. 34 

0.35 

0. 62*2 

q(D) 

-1.03 


-0. 12 

-0. 24 

-0.22 

-0. 35*2 

q < BaO ) 

0.00 

— 

♦ 0. 06 

♦1.27 

♦ 1.25 

♦0.87*2 

q<V» ) 


-0.21 

-1.05 

♦ 0. 06 

-0. 07 

m 

o 

o 

q ( W i i ) 

— 

-0. 30 

♦ 0. 19 

-0. 33 

-0.23 

-0.09 

qIU.,,1 


-0.53 


-1.00 

-0.87 

-0.86* 

q(V w) 

.... 




-0.08 


ft (Debye) 

13.8 

25. 1 

14.9 

37.5 

38.2 

24.6*2 

E dipole ( e V ) 

5.5 

3. 1 

4. 4 

6. 4 

5.0 

5.1*2 

E irn fh ( eV ) 

... 

2.5 

(2.2) 

5.2 

4.9 

3.4*2 

* -0. 16 for 
-0.38 for 

U * (center ) , +0. 

U i » i ( cent e r ) , - 0. 

19 for V, 
48 for W, 

(periphery), 
in (perl phery ) . 



suits in very large dipole moments of about 38 Debyes for 
the larger single-BaO-adsorbate clusters. The dipole mo- 
ment for BaO/W 9 is 1.5 times larger than for Ba/W 9 , 
assuming the same Ba-surface distance. The presence of 
oxygen thus strongly enhances the charge transfer from 
Ba into the W substrate, which leads to the observed low- 
ering of the work function [22], [23]. The large initial 
dipole moment for BaO/W ( 100) of about 38 Debyes per 
BaO adsorbate will be significantly reduced in a com- 
pleted overlayer because of repulsive dipole-dipole inter- 
actions. The result for Ba 2 0 2 /W 14 with a net charge of 
+ 1.2 electrons on Ba and an effective dipole of 25 De- 
byes per BaO adsorbate shows the expected strong depo- 
larization of interacting dipoles. 

The maximum coverage on W( 100) and real cathode 
surfaces is probably 1 /2 or less, and the BaO-BaO sep- 
aration is therefore likely to be larger than in the W, 4 clus- 
ter. On the other hand, there will be more nearest neigh- 
bors that may offset the possibly too large depolarization 
resulting from one neighbor at a closer distance. The net 
charge on Ba and the dipole moment for the single-BaO- 
adsorbate clusters with a formal coverage 1 /4 are there- 
fore expected to be reduced for larger coverages on real, 
polycrystalline cathode surfaces, and the results for 
Ba 2 0 2 /W 14 are certainly considered to be more realistic. 
Analyses of Auger spectra suggest that Ba is in a mixed 
metal-metal oxide valence state [24]. 

An additional test for the accuracy of the calculated 
charge distribution in the surface clusters is provided by 
Ba binding energies obtained from simple electrostatic and 
image charge concepts, which are also given in Table II. 
The dipole and image charge binding energies are defined 
as, respectively, 


^dipole ^Ba /^dipole *?Ba / M 


and 


^im.ch. 




Ba-surface* 


Atomic units are used with q in electrons, R in Bohrs and 
H in a.u. (1 A,%/0.5292 Bohrs, 1 Debye /^'2. 542 a.u,, 
and 1 eV/%^27.21 Hartrees); /? Ba surface = 6.0 Bohrs. 

The validity of the dipole expression for the binding 
energy in ionic systems has been tested for molecular BaO 
at the equilibrium separation, R e = 1.94 A [25]. Assum- 
ing point charges, the experimental dipole moment of 7.93 
Debyes corresponds to a transfer of 0.85 electrons from 
Ba to 0. The resulting dipole binding energy is 5.4 eV, 
which is in reasonable agreement with the experimental 
value for molecular dissociation of 5.8 eV. Relativistic 
Xa calculations give binding energies of 6. 1 eV for R e and 
5.5 eV for the stretched molecule with R = 4.9 Bohrs. 

The calculated Ba binding energy for the larger single- 
BaO-adsorbate clusters is of the order of 6 eV and reduced 
by dipole-dipole interaction to 5. 1 eV in Ba 2 0 2 /W 14 . The 
latter value is in excellent agreement with the experimen- 
tal value of 4.8 ± 0.2 eV for Ba desorption [22]. If simple 
image charges for q Ba are assumed inside the substrate and 
Rimxh. = 2 ^Ba-surface = 12 Bohrs are used (rather than the 
actual center of the negative charge via the calculated di- 
pole moment), the Ba binding energies are reduced by 
about 1 eV. The reduction is related to the fact that the 
actual dipole length R dipo i e = fx jq Bd is smaller than R im ch 
because of the presence of a small negative charge on ox- 
ygen and screening effects inside the W substrate. The 
calculated dipole binding energy for Ba/W 9 at R = 6 
Bohrs is 3. 1 eV, which may be compared with the exper- 
imental value of 3.8 eV for Ba desorption from a clean 
tungsten surface [22]. 

In summary, it is found that the electronic structure of 
the BaO/W emissive layer appears to be reasonably well 
described by the present calculation. The net positive 
charge of 1.2 electrons on Ba for the largest cluster gives 
rise to a Ba binding energy of about 5 eV in agreement 
with experiment. The net negative charge on O is less 
than 0.5 electrons. There are, however, indications from 
the comparison of calculated and experimental densities 
of states for BaO/W (100) that the O-W distance may be 
larger than assumed. The comparison of the Ba/W and 
BaO/W interaction shows that the O/W electronic struc- 
ture provides for a substantially increased charge transfer 
from Ba into the metal and a larger binding energy. The 
resulting very large surface dipole is responsible for the 
lowering of the work function and related improvement 
in the emission properties of the cathode. The strong ionic 
bond in the Ba-O/W complex explains the high thermal 
stability of the surface. 


V. BaO/W-Os Model Results 

As a first step toward the theoretical investigation of 
alloy cathode surfaces, different bcc unit-cell substrate 
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clusters have been considered in which W atoms were se- 
lectively replaced by Os. All other computational param- 
eters and the adsorption geometry are kept the same in 
order to emphasize the effects due to the atomic substi- 
tution and eliminate geometrical and cluster size and shape 
effects. 

Three different BaO/W-Os clusters are derived from 
W 9 by replacing either the W b W H , or W m atoms by Os, 
which leads to substrate clusters denoted as Os 4 WjW 4 , 
W 4 0s!W4, or W 4 W,Os 4 . In the W 2 0s 2 W 1 0s 2 W 2 cluster, 
two first- and two third-layer W atoms in the [ 110] and 
[ 110] directions, respectively, are replaced by Os. 

The calculated charge distributions, dipole moments, 
and Ba binding energies for these clusters are given in 
Table III together with the results for BaO/ W 9 and the Ba 
bridge-bonded Ba0 2 /Os 9 cluster. The comparison of the 
data for the mixed BaO/W-Os clusters with BaO/W 9 
shows two significant results: First, the net charges on Ba 
and O are very similar, with a maximum variation of 0. 1 
electrons only. Second, very large charge redistributions, 
however, occur in the substrate which affect in a system- 
atic way the cluster dipole moments and Ba dipole bind- 
ing energies. This situation is illustrated in Fig. 4. 

The charge redistribution in the substrate has to be at- 
tributed to the electronegativity difference between W 
(1.40) and Os (1.52) [26]. The higher electronegativity 
value for Os indicates that Os has a larger tendency to 
attract electrons of chemical bonds than W, which results 
in a net charge transfer from W to Os. This has the fol- 
lowing consequences for the surface dipoles as qualita- 
tively indicated with — / + charges and dipole arrows at 
the bottom of Fig. 4: If all Wj atoms are replaced by Os, 
charge is transferred from W n and W m to Os] opposing 
the Ba-substrate charge transfer, and the BaO/W 9 (ref- 
erence) dipole is lowered by 8 Debyes. If all W IH atoms 
are replaced by Os, charge is transferred from W t and W n 
to Os [H in the direction of the Ba-substrate charge trans- 
fer, and the BaO/W 9 dipole is increased by 4 Debyes. If 
W n or half of the W atoms in the first and third layers are 
replaced by Os, the resulting cluster dipoles are within 
about 1 Debye of BaO/W 9 . 

If Os is present in the first layer, as in 
BaO/ W 2 Os 2 W , Os 2 W 2 or BaO/Os^W^, the net charge 
on barium is increased relative to BaO/W 9 , which results 
in a larger Ba dipole binding energy. If Os is not present 
in the first layer, the Ba binding energy is smaller than in 
BaO/W 9 . In the hep Ba0 2 / Os 9 cluster, where Ba is closer 
to the surface, the net positive charge on Ba is reduced 
from about 1.5 to 1.3 electrons and the dipole moment is 
reduced to almost one half of the BaO/W 9 value, which 
results in a net increase in the Ba dipole binding energy. 

The present theoretical results for simple model alloy 
systems suggest that the following general trends should 
hold: 1) For similar BaO geometries, the presence of Os 
in the first layer leads to a smaller surface dipole and larger 
Ba binding energy compared to pure W. 2) If Ba is bridge 
bonded to O and consequently closer to the surface than 



y. CD) 36 37 30 37 k2 

E dip (ew 7,0 6,11 9,2 6J 5,6 

Fig. 4. Charge distribution in BaO/(W-Os) 9 clusters. The +/- signs in- 
dicate the charge redistribution in the substrate layers and the arrows the 
dipole contribution relative to the BaO/W 9 reference cluster. 

TABLE III 

Net Charges ?< Electrons), Dipole Moments y , and Ba Dipole 
Binding Energies E dipidt . for BaO/W-Os Clusters 
(R(Ba-O) — 4.9 Bohrs ~ 2.6 A.) 


Property 






B»0 a / 
Os«( hep) 

U 3 0s 3 U,0s a U 

2 u* 

Os«U,V« 

U*0s i V« 

VUU ,Ds 4 

q(Ba) 

>1.54 

>1.51 

>1.58 

>1.48 

♦1 . 50 

>1.31 


- 

I 





q<0) 

-0.22 

-0.24 

-0.17 

-0.21 

-0. 24 

-0. 36 « 2 

q(B^D) 

* >1.32 

>1.27 

>1.41 

>1.27 

>1.26 


q" f W i ) 

>0. 38 

>0.06 


>0. 15 

>0. 29 


q ( 0 * , > 

-0.63 


-0.64 



>0, 05 

q(U,>) 

-0. 11 

-0. 33 

-0. 09 


-0.09 






-0. 54 


-0. 15 

q(U in ) 

-0.04 

-1.00 

-0.48 

-0.08 



q ( Os z i x ) 

-O. 92 




-1.46 

-0. 25 

fl(Ol.u) 






-0.24 


y (Debye) 

36.2 

37.5 

29. 5 

36.6 

41.8 

19.7 

^dipole I*V) 

7.0 

6.4 

9.2 

6. 1 

5.6 

7.9 


in the on-top position the surface dipole decreases. 3) For 
different alloy substrates, the largest changes in the charge 
distribution occur in the substrate, which can either com- 
pete with or enhance the Ba-substrate charge transfer, de- 
pending on the location of the different metal atoms in the 
first few layers. 

A higher Ba binding energy is experimentally observed 
for osmium-ruthenium coated W surfaces as present in 
commercial M-type thermionic cathodes. The binding en- 
ergy increases from 4.8 eV for B-type to 5.4 eV for M- 
type cathodes [27]. The difference correlates very well 
with the calculated increase from 6.4 eV for BaO/W 9 to 
7.0 eV for BaO/W 2 Os 2 W 1 Os 2 W 2 at presumably lower 
coverages. The calculated Ba binding energy on a pure 
Os substrate in Ba0 2 /0s 9 is 7.9 eV. 

The question of the lowering of the work function in 
alloy cathodes cannot easily be addressed from the pres- 
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ent results for individual surface dipoles because the work 
function is an extended surface property. The minimum 
work function results from the competition of an increase 
in the total surface dipole due to additional adsorbate di- 
poles and the depolarization that occurs with increasing 
coverage. Simple models for work function versus cov- 
erage curves suggest that the maximum work function 
change is determined by the magnitude of the initial sur- 
face dipole, while the optimum coverage is determined by 
the strength of the depolarization [28]. A larger surface 
dipole would therefore be indicative of a lower work 
function, and the on-top arrangement should theoretically 
lead to a lower work function at smaller coverages. Be- 
cause, however, the work function versus coverage curve 
for the BaO/W system does not exhibit a minimum but 
continues to decrease with increasing coverage [23], the 
bridge-bonded BaO overlayer with smaller individual di- 
poles could lead to a lower work function at higher cov- 
erages. 

The difference between B- and M-type cathodes has 
been related by Norman et ah [6], in part, to the different 
on-top (BaO/W) versus bridge-bonded (BaO/W-Os) 
geometries, which were questioned, however, by Shih et 
al. [7] and are inconsistent with the model for the bonding 
proposed by Green [29]. Drawing on theories of hetero- 
geneous catalysis, Green had postulated that the Ba-0 
surface bond is tilted on the tungsten substrate and linear 
on alloy substrates. From the present model results for 
linear geometries, the lower work function for M-type 
cathodes could possibly be explained by a higher final 
coverage due to reduced dipole-dipole interaction and/or 
higher Ba binding energy. 

VI. Summary and Conclusions 

Detailed analyses have been presented for the elec- 
tronic structure of the BaO/W surface complex as ob- 
tained from self-consistent relativistic scattered-wave 
cluster calculations. The effects of mixed W-Os alloy sub- 
strates are investigated with computational methods for 
the first time. Calculated properties include electronic 
energies, densities of states, electron densities, charge 
distributions, surface dipoles, and barium binding ener- 
gies on different model substrates. 

Barium is found to be heavily oxidized in the BaO/W 
surface complex, with an oxidation level greater than 1 , 
while both oxygen and the tungsten substrate are in a par- 
tially reduced chemical state. The net negative charge on 
oxygen is less than 0.5 electrons and thus surprisingly 
small. The remaining negative charge balancing the large 
positive charge on barium is found in the form of a 
screened image charge inside the substrate cluster. The 
bonding is largely ionic between Ba and the O/W sub- 
strate, while oxygen forms covalent bonds with tungsten. 
The chemical structure of the surface is therefore char- 
acterized as a Ba salt of the O/W substrate. These theo- 
retical findings are consistent with available experimental 
data which describe the barium oxidation state as inter- 
mediate between metal and oxide, oxygen in an environ- 


ment different from bulk BaO, and the W substrate as not 
being oxidized [29]. 

The presence of oxygen leads to a significant enhance- 
ment of the surface dipole relative to Ba/W. The theo- 
retical population analyses indicate that the O 2 p orbitals 
are strongly mixed with W 5 d orbitals and thus provide 
for an effective charge transfer from Ba into the W sub- 
strate. The large charge transfer leads to a very large ini- 
tial surface dipole of up to 40 Debyes for the assumed 
geometry of Ba on top of O in the fourfold hollow site of 
W( 100). The repulsion between interacting adsorbate-in- 
duced dipoles, however, causes substantial depolariza- 
tion. The analysis of the charge distribution shows that 
the relevant surface dipole is the Ba-substrate and not the 
Ba-0 dipole. The calculated Ba binding energy for an in- 
dividual BaO adsorbate is about 6 eV. For strongly inter- 
acting BaO adsorbates this value is reduced to 5 eV, which 
is in excellent agreement with the experimental energy for 
Ba desorption from BaO/W surfaces. 

The comparison of theoretical and experimental O 2 p 
binding energies in the BaO/W ( 100) system suggests that 
the oxygen-surface distance of 0.6 A as derived from a 
fit of surface models to experimental SEXAFS data [17] 
may be too small. For a larger oxygen-surface distance in 
BaO/W, the net charge transfer from W to O increases 
and the surface dipole decreases [21]. The surface dipole 
also decreases when barium is bridge-bonded to two ox- 
ygen atoms (and is consequently closer to the surface) as 
shown with the model for adsorption on a hexagonal- 
close-packed osmium cluster. 

The analysis of the charge distributions in BaO/W-Os 
model alloy clusters with identical adsorbate and substrate 
geometYies shows similar net charges on Ba and O but 
very large charge redistributions inside the substrate clus- 
ters. The charge redistribution is attributed to the higher 
electronegativity of Os relative to W and varies system- 
atically with the positions of the Os and W atoms in the 
clusters. The surface dipole may either increase, de- 
crease, or remain similar to BaO on pure W. If Os is pres- 
ent in the surface layer, the charge transfer from Ba to the 
substrate and the Ba binding energy increase relative to 
W. The calculated Ba binding energy for a more or less 
random W-Os substrate cluster is found to be 0.6 eV 
higher than for the equivalent pure W cluster, a result that 
is in striking agreement with the experimentally observed 
increase in the energy for Ba desorption from B- to M- 
type surfaces. 

The results from the present theoretical investigation 
suggest that for assumed similar BaO geometries the ob- 
served lowering of the work function and related emission 
enhancement for alloy versus tungsten cathode surfaces 
could be explained by higher Ba coverages due to higher 
Ba binding energies. The higher Ba binding energies re- 
sult from an increased charge transfer from Ba to the sur- 
face and the charge redistribution within the first substrate 
layers due to electronegativity differences between the al- 
loy component metals. 

On the basis of the electronegativity argument for Os 
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and W, one might speculate that a Pt-W surface should 
also show an improved emission relative to W, since the 
electronegativity for Pt (1.44) is also slightly higher than 
for W (1.40) [26]. This is, of course, well known to be 
not the case. A significant difference between Pt and the 
emission enhancing transition elements Re, Os, Ru, and 
Ir is the fact that Pt has an almost completely filled 5 d 
shell and forms intermetallic compounds with Ba [30]. 
The electronegativity argument, therefore, is not equally 
applicable to Pt, and these other properties are obviously 
much more important. 

A consideration of different bonding geometries in the 
Ba-0 overlayer shows that the surface dipole for an in- 
dividual adsorbate decreases when changing from on-top 
to bridge-bonded BaO (with Ba closer to the surface). 
The larger dipole for the on-top orientation would favor a 
lower work function at smaller coverages, if such an 
overlayer actually forms and Ba and O are not initially 
adsorbed in different sites. Because of the strong repulsive 
interaction between adsorbate-induced dipoles, however, 
the bridge-bonded BaO overlayer with smaller individual 
dipoles would tend to be more stable at higher coverages 
and could thereby lead to a lower effective work function. 
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